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Introduction 


-^The  work  of  this  contract  was  intended  to  develop  the  atmospheric 
model  which  integrates  the  emission  and  the  scattering  term  in  a  combined 
form  and  to  provide  analysis  of  the  earth  atmospheric  emission  data  gathered 
by  the  SPIRE  payload  in  these  two  terms. 

Our  work  performed  throughout  the  entire  period  can  be  divided  into 
five  areas: 

(1)  incorporation  of  various  computer  programs  designed  to 
calculate  the  atmospheric  infrared  emission  level  under  a 
specific  condition; 

(2)  improvement  of  the  algorithm  employed  in  these  programs; 

(3)  development  of  an  atmospheric  infrared  radiance  model  which 
unifies  the  emission  and  scattering  term; 

(U)  application  and  evaluation  of  the  developed  model  to  analysis 
of  the  SPIRE  data;  and 

(5)  evaluation  of  some  preliminary  test  data  for  the  coming  AFGL 

balloon-borne  atmospheric  emission  measurement  experiment. 

The  SPIRE  experiment^ was  carried  out  during  1978,  and  the  collected 

data  were  made  available  to  us  in  late  1978.  Our  effort  during  the  early 

phase  of  the  contract  period,  i.e.,  prior  to  the  SPIRE  launching,  was 

focused  in  the  area  of  (l)  and  (2).  Our  Scientific  Report  No.  1  summarizes 
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our  progress  made  during  the  early  phase.  The  present  report  will  describe 
our  work  covered  in  the  entire  period,  with  a  special  emphasis  on  the 
progress  made  since  the  last  report. 

>n 

Upon  inspection  of  the  SPIRE  data,  it  became  obvious  to  us  that  the  _ 

existing  computation  programs  are  much  more  elaborate  than  required  in  the  — 
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analysis.  A  quick  check  on  the  observed  scattering  radiance  level  revealed 
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that  the  single-scattering  theory  with  the  Rayleigh  particles  produced  a 

result  adequate  for  the  analysis.  The  Monte  Carlo  code  developed  by 
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Blattner  and  his  co-workers  is  found  unnecessary  for  the  analysis,  because 
the  computation  based  on  the  single  scattering  process  produced  a  result 
which  agrees,  reasonably  well,  with  the  SPIRE  observed  data.  The  existing 
computation  programs  treat  the  scattering  process  entirely  separately  from 
the  absorption-emission  process.  We  find  them  rather  inadequate  for 
explaining  why  some  molecular  bands  change  their  absorptive  behavior 
against  the  scattering  background  to  become  emissive  as  the  altitude 
increases.  A  theory  which  unifies  both  processes  in  the  radiative 
formulation  is  definitely  needed  for  the  analysis.  A  limit  placed  on  time 
restricted  us  to  develop  a  computation  algorithm  which  generates  the 
spectral  radiance  level  as  a  function  of  frequency  and  an  observation 
condition  in  an  automatic  fashion.  We  feel  that  we  are  on  the  right  track 
for  seeking  a  solution  to  the  problem.  Nonetheless,  we  came  to  understand 
the  basic  radiative  mechanism  of  the  upper  atmosphere  with  an  extent  that 
we  can  place  a  proper  perspective  to  the  analysis  of  the  SPIRE  data. 

Radiative  Process 

The  radiative  transfer  in  the  atmosphere  is  traditionally  analyzed  in 
treating  the  infrared  radiation  field  as  a  continuous  flow  of  radiative 
energy.  The  interactions  between  the  radiation  field  and  the  atmospheric 
gas  is  formulated  in  an  integral  equation.  Several  computer  programs  have 
been  written  for  analysis  of  the  radiation  transfer  problem.  Our  approach 
taken  in  the  present  study  is  slightly  different  from  the  traditional 
analysis.  We  tried  to  analyze  the  radiative  transfer  in  the  upper  atmosphere 
in  following  the  interactive  process  between  the  photons  and  molecules. 


The  photons  incident  into  a  layer  of  atmosphere  would  go  through  one 
of  three  processes:  (l)  they  may  be  absorbed  by  the  atmospheric  molecules 
or  atoms;  (2)  scattered  by  the  Rayleigh  particles  or  by  the  aerosols;  and 
(3)  left  unperturbed.  We  can  formulate  these  interactions  between  the 
photon  field  and  the  molecular  system  in  terms  of  the  density  of  these  two 
interacting  partners  and  a  quantity  which  indicates  a  strength  of  the 
interaction. 

The  interaction  between  the  photons  and  the  molecules  (or  the  atoms) 

8 

in  the  absorption  process  is  specified  by  Einstein's  formulation.  The 

number  of  molecular  transitions  in  unit  time  per  unit  volume,  when 

interacting  with  the  photons  of  a  frequency  v,  is  given  by 

N,  =  (N  B  -N  B  )p( v) .  (1 

t  £  £u  u  u£ 

In  this  expression,  the  following  notations  are  used: 

N^:  number  of  molecules  per  unit  volume  in  lower  state  £ 

N^:  number  of  molecules  per  unit  volume  in  upper  state  u 

B  :  Einstein  coefficient  for  induced  absorption 

B  :  Einstein  coefficient  for  induced  emission 

ujl 

p(v):  spectral  density  of  the  radiation  field  in  unit  volume  at 
frequency  \> 

The  density  p(v)  is  related  to  the  photon  density  Np  given  in  unit  volume 
per  second  by 

fp(v)  =  Nphv,  (2 

where  h  is  the  Planck  constant  and  c  is  the  speed  of  light.  Combining 


Eqs.  (l)  and  (2),  the  molecular  transitions  Np  is  given  by 


IT  =  ( N.B  -N  B  J 
t  £  Jlu  u  u£ 


UN  hv 


This  expression  can  be  rewritten  by  using  the  molecular  absorption 


strength  S  (cross-section) 
B  £u 


where 


r  =  1*N  N  S„  , 
t  p  o  lu 


U) 


0  hv(N  B  -N  B  „)  .  N 

S„  =  2,  2u  u  u«,  .  hv,,  , ,  uN 

£u  - - -  =  B.  (l-TT-J 

N  c  c  £u  N 

o  l 


(5) 


with  N  the  number  of  molecules  in  unit  volume.  In  deriving  Eq.  (5),  it  is 
o 


as  siomed  that 


N  =  N  . 
o  & 


(6) 


The  Eq.  ( U )  derived  above  specifies  the  number  of  molecular  transitions 
in  the  photon  density  N  ,  the  molecular  density  Nq,  and  the  molecular 
absorption  strength  S  which  is  experimentally  determined  in  the  laboratory. 

The  scattering  process  can  be  formulated  in  the  same  way.  A  traditional 
approach  is  somewhat  different  in  that  the  incident  radiation  field  is  given 
in  the  irradiance;  the  power  scattered  toward  the  direction  4>  by  a  molecule 
is  given  by 


0  ( 4> )E  =  1(b), 
m  o 


(T) 


where  E  is  the  irradiance  of  the  incident  plane  wave.  We  can  formulate  the 
o 


incident  irradiance  by 


E  =  N  hv, 
o  p 


(8) 


again  using  N  ,  the  photon  density  in  unit  volume  per  unit  second.  Now  we 
find  that  the  power  scattered  by  a  molecule  is  given  by 


1(b)  =  N  hv  a  (b). 
p  m 


(9) 


The  radiative  power  scattered  toward  the  direction  b  per  unit  solid  angle  by 
an  ensemble  of  molecules  N  in  unit  volume  is  given  by 


I'(b)  =  N  hv  N  a  (b). 
o  pm 


(10) 


Thus  we  get  the  desired  form  for  the  number  of  photons  scattered  toward  the 


-  — S+  VS 


direction 


N  (<j>)  =  N  N  o  U). 
p  o  p  m 

For  the  Rayleigh  scattering,  the  scattering  cross-section  is  given  by 

2 


(11) 


o  <♦)  -**£( 

m 


n2-l 


N2c4  n2+2 


)  sin2 


(12) 


Utt2  V2 


(n-l)2  sin2 


N2c2 

It 

where  N  is  the  number  of  molecules  in  a  unit  volume  at  S.T.P. 

The  energy  absorbed  by  a  molecule  from  the  radiative  field  is  sooner 
or  later  released  from  it  in  a  de-excitation  process.  In  the  upper 
atmosphere,  the  de-excitation  is  done  by  the  radiative  process.  The 
photon  absorbed  by  a  molecule  is  re-emitted  to  the  radiation  field 
isotropically.  If  the  original  absorption  occurs  between  the  ground  state 
and  the  first  excited  vibrational  state,  there  is  no  significant  change  in 
the  photon  energy  of  the  absorbed  and  of  the  re-emitted.  In  the  lower 
atmosphere,  a  local  thermal  equilibrium  is  established;  the  energy  absorbed 
by  a  molecule  is  dissipated  into  the  kinetics  energy  shared  in  the  molecular 
ensemble.  The  interaction  between  the  molecular  ensemble  and  the  radiation 
field  under  the  J coal  thermal  equilibrium  condition  is  treated  below  in 
calculation  of  the  photon  density.  The  de-excitation  under  this  condition 
is  done  predominantly  by  the  collisional  process.  In  a  transient  region 
from  the  lower  to  the  upper  atmosphere,  the  molecular  system  de-excites  by 
both  processes. 


Photon  Density  In  Atmosphere 

There  a'-e  two  radiation  sources  which  must  be  considered  in  the  infrared 
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radiative  process  of  the  atmosphere:  the  solar  radiation  penetrating  the 
atmosphere  from  the  top,  and  the  radiation  going  upward  to  escape  to  the 
outerspace  froin  the  earth's  surface.  The  energy  density  p(v)dv  of  the 
black-bodv  at  temperature  T  is  given  by 


o ( v ) dv  = 


Snhv^  dv 


a  hv 
kT 

e  -1 


(13) 


Hie  photon  density  in  unit  volume  per  unit  second  per  wavenumber  interval 
da  emitted  by  the  same  blackbody  radiator  is 

r2 

K  da  =  r-r—  p(v)dv 
p  4hv 


2ttv^  da 


2na2c  da 


hoc 

kT  , 
e  -1 


(1U) 


hca 

kT 


-1 


For  the  telluric  radiation,  the  photon  density  is  calculated  by  taking  T  at 
a  value  between  225°K  and  300°K.  For  the  solar  radiation,  it  is  calculated 
by  the  same  expression  with  an  attenuation  factor; 


N  =  (— ) 

P.  V 

kT  , 
e  -1 


(15) 


where  r  is  the  solar  radius  and  R  is  the  solar  distance,  T  is  about  6000°K. 
s  s 

The  front  factor  (rs/R  )  is  approximately  2.2  x  10  .  In  Figure  1,  the 
photon  density  for  three  cases  is  plotted  as  a  function  of  the  wavenumber  a: 
(l)  the  solar  radiation  calculated  by  Eq.  { 1 5 ) ,  (2)  a  blackbody  at  T  =  300°K, 
and  (3)  a  blackbody  at  T  =  225°K.  The  telluric  photon  density  in  the  upper 
atmosphere  is  affected  by  the  optical  thickness  of  the  lower  atmosphere.  In 
some  spectral  ranges,  there  is  no  atmospheric  molecular  absorption,  and  the 
photon  density  for  the  300°K  blackbody  remains  unaffected  from  the  earth's 
surface  to  the  top  of  atmosphere.  In  a  region  where  the  atmospheric 
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absorption  is  extreme,  the  photon  density  in  the  upper  atmosphere  is 
controlled  by  the  energy  density  of  the  blackbody  radiation  at  the  tropopause 
temperature  T  =  225°K.  We  can  observe  in  Figure  1  that  the  molecules  and  the 
atoms  in  the  upper  atmosphere  are  excited  by  two  types  of  radiations  which 
have  a  vastly  different  spectral  characteristic.  The  telluric  photon  density 
remains  unaffected  by  the  solar  daily  motion.  The  photons  in  a  region  below 
2000  cm  are  predominantly  telluric,  while  those  in  a  region  above  2000  cm  ^ 
in  the  daytime  are  solar. 

,  3-0 

The  emission  due  to  the  CO^  transitions,  one  for  the  (01101-00001)  band 

at  670  cm  (the  band)  and  another  for  the  ( 00011-00001)  band  at  2300  cm 

(the  band),  is  ready  for  a  qualitative  analysis  using  the  photon  density 

calculated  for  the  curves  shown  in  Figure  1.  (A  more  detailed  quantitative 

analysis  will  follow  later  in  the  report. )  Once  accepting  the  photon 

spectrum  shown  in  Figure  1,  we  find  that  the  photons  available  for  exciting 

the  band  are  telluric.  The  temporal  characteristics  of  the  telluric 

photons  remain  independent  of  the  daily  cycle.  Thus  the  ( 01101-00001 ) 

band  emission  at  670  cm  1  must  remain  unchanged  between  the  day  and  the 

night  time.  The  SPIRE  data  shown  in  Figure  2  for  this  CO^  band  confirms  a 

stationary  nature  of  the  emission.  The  00^  ( 00011-00001 )  band  at  day 

side  is  excited  by  the  solar  photons,  and  by  the  telluric  photons  at  night 

side.  Since  the  transition  is  very  strong,  the  telluric  photons 

available  in  the  upper  atmosphere  are  those  of  a  225°K  blackbody.  The 

difference  in  these  two  densities  between  the  solar  and  the  telluric  is  two 

orders  of  magnitude.  Thus  the  radiance  level  of  the  transition  should  undergo  a 

change  of  about  2  orders  of  magnitude  between  the  day  and  the  night.  The 

SPIRE  data  in  Figure  3  indicates  that  the  band  emission  primarily  under¬ 
goes  the  diurnal  change  described  above.  Upon  a  careful  inspection  of  the  data, 
the  radiance  level  at  nighttime  is  brighter  than  what  can  be  predicted  by  the 
telluric  photon  level. We  will  make  an  analysis  on  this  anomaly  later  in  the 
report. 
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Figure  3 
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Model  for  Calculating  Atmospheric  Radiance  Level 

The  expressions  obtained  above  in  Eqs.  (1*)  and  (ll)  are  concerned  with 
an  individual  interacting  process  in  which  a  single  photon  is  either  absorbed 
or  scattered  by  a  molecule.  In  this  section  we  will  develop  a  model  which 
relates  the  atmospheric  radiance  level  to  the  basic  interaction  process 
given  by  these  expressions. 

The  energy  absorbed  by  a  molecule  in  interaction  with  the  photon  field 
must  be  released  in  a  de-excitation  process.  In  the  upper  atmosphere,  the 
only  process  available  for  de-excitation  is  a  radiative  process.  A  question 
which  must  be  raised  and  answered  for  such  a  radiative  excitation-de¬ 
excitation  process  concerns  an  establishment  of  equilibrium  population 
density  among  various  molecular  states.  Degges  and  other  workers  developed 

the  computation  program  capable  of  giving  an  estimate  of  the  radiative  level 

11,12 

in  the  normal  atmospheric  environment  at  high  altitude.  We  are  truly 
thankful  for  the  existence  of  these  programs  which  provided  a  valuable  help 
to  develop  our  insight  into  the  problem. 

We  modified  the  Degges  program  to  suit  our  needs,  and  to  improve  the 
computation  efficiency  in  our  CYBER  system.  A  brief  description  of  our 
effort  along  this  direction  will  be  given  later  in  this  report,  together 
with  the  listing  of  our  modifications.  The  atmospheric  model  is  judged  on 
its  ability  for  producing  a  reasonable  estimate  of  the  infrared  radiance 
level;  in  particular,  it  is  questioned  on  precision  and  practicality. 

The  computation  scheme  to  obtain  the  radiance  level  must  find  the 
equilibrium  population  distribution  which  is  established  among  the 
vibrational  states  when  the  molecular  system  is  subjected  to  balance  with 
the  photon  field.  The  precision  needed  in  the  computation  is  somewhat 
controlled  by  a  spectral  resolution  of  the  result.  The  emission  involving 
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the  hot-band  transition  can  be  distinguished  only  with  an  improved  spectral 
resolution,  because  it  is  usually  buried  under  the  fundamental  transition. 

We  decided  to  ignore  these  hot-band  transitions  from  our  model  in  a  belief 
their  contribution  is  not  significant  to  the  computation.  The  assumptions  made 
in  our  model*  are  certainly  open  for  debate  and  all  criticisms  to  our  view¬ 
point  are  welcome.  Our  view  that  all  transitions  observed  in  the  data 
occurred  from  the  ground  state,  in  turn,  simplified  the  calculation  to  a 
great  extent.  A  direct  correlation  was  resulted  in  our  formulation  between 
the  radiance  level  for  a  particular  tangent  height  and  the  molecular  density 
of  the  concerned  molecular  system  at  this  altitude. 

The  total  energy  emitted  by  the  molecules  in  a  unit  volume  through  the 
radiative  de-excitation  is  equal  to  the  energy  absorbed  by  the  same 
molecules  if  that  is  the  only  process  feasible.  The  absorbed  energy  is 
given  by  the  number  of  transition  times  the  photon  energy:  b'  =  hvN^. 

The  energy  released  in  the  radiative  de-excitation  is  emitted  isotropically. 
Thus  the  energy  re-emitted  per  a  unit  area  over  a  unit  solid  angle  is  given 
by:  b  =  b'Air  =  hvN^/liTr.  The  photons  released  from  a  molecule  are  very 
unlikely  absorbed  again,  since  the  atmosphere  is  optically  thin.+  Even  if 
they  are  absorbed,  they  are  re-emitted  in  the  photon  field.  The  photons 
which  reach  a  sensor  of  a  unit  area  with  a  field  view  of  a  unit  solid  angle 
are  computed  simply  by  accumulating  the  contribution  from  each  unit  volume 
which  falls  into  a  line  of  sight.  Since  the  optical  etendu  preserves,  the 
summation  which  we  must  take  in  the  calculation  is  equal  to  that  of  each 
contribution  b  taken  over  a  column  of  a  unit  cross-section  extending  from 
the  sensor  to  the  space  all  the  way  along  the  line  of  sight.  The  radiance 
level  per  unit  solid  angle  integrated  over  the  band  is  then  given  by 

*See  the  Appendix  D. 

tThe  absorption  along  the  limb-view  path  of  100  km  tangent  height  is  less 
than  5$  for  the  ( 00011-00001 )  band  of  C02. 


where  N  is  the  total  number  of  molecules  counted  along  the  line-of-sight 


per  unit  area; 


N  =  /N  di  . 
o 


(17) 


The  density  of  atmospheric  molecules  falls  off  approximately  exponentially 
with  the  altitude.  If  we  can  assume  the  density  N(h)  as  a  function  of 
altitude  h  given  by  ^  ^ 


JV(h)  =  N(h  )  e 
o 


(18) 


with  a  characteristic  constant  H,  the  total  number  of  molecules  between  two 


altitudes  h  and  h  is  given  by 

o 


h-h 


N(h,h  )  =  fh  N( h)  dh  =  N{h  )  fh  e  H  dh 
o  ,  o 

h  h 

o  o 


(19) 


=  N(h  )  H 
o 

For  the  limb  view  calculation,  we  need  to  calculate  a  distance  l 
between  two  altitudes,  the  tangent  height  x  and  a  height  above  by  H.  The 
distance  l  is  calculated  by 

l  =  (x+H+p)  Jl  -  (^J^?  =  ^SH)  (x+p)  (20) 

where  p  is  the  earth  radius  (-  6360  km).  To  calculate  the  total  density  N 
along  a  line  of  sight  to  the  tangent  height  x,  we  use  the  approximation 
Eqs. (19)  and  (20) : 

N  =  N  (x)  2  /2(H)  (x+p)  (21) 

o 

where  N^(x)  is  the  molecular  density  at  altitude  x.  With  the  approximate 
total  density  N  derived  in  Eq.  (21)  the  limb-view  radiance  at  the  tangent 


height  x,  is  given  by  an  expression 

B  =  —  S  „  N  (x)  2  /2H(x+p )  N  .  (22) 

it  uZ  o  p 

The  expression  given  by  the  equation  above  requires  an  estimate  of  the 
photon  density  N^.  The  value  shown  in  Figure  1  is  given  in  a  unit  wavenumber 
interval.  The  photons  which  contribute  to  the  transition  are  those  within 
the  width  of  the  vibrational-rotational  lines.  If  we  take  the  photon  density 
to  be  the  value  over  the  band,  say  100  cm  ^  for  the  CO^  ( 01101-00001 ) 
transition,  we  obviously  get  an  erroneous  result.  The  effective  photon 
density  must  be  adjusted  by  multiplying  an  effective  number  of  the  rotational 
lines  within  the  band  times  the  line  width,  which  we  can  safely  assume  the 
Doppler  width  for  those  in  the  upper  atmosphere. 

C0o  Bands,  ( 01101-00001 )  and  ( 00011-00001 ) 

The  CO^  molecule  is  relatively  stable  in  the  atmosphere,  because  neither 
electronic  transitions  nor  photo-dissociation  occurs  in  the  entire  spectral 
range  below  the  Schuman-Runge  Continuum  of  0^-  The  mixing  ratio  remains 
constant  in  the  atmosphere  below  120  km  where  the  0^  absorption  is  very 
effective  in  blocking  those  high  energy  photons  in  the  extreme  uv  region. 

As  a  result,  the  equilibrium  population  distribution  is  established 
among  various  vibrational  levels  in  the  electronic  ground  state  with 
insignificant  perturbation  from  higher  excited  electronic  states.  The 
equilibrium  is  achieved  in  balance  with  the  infrared  photon  field.  Under 
this  condition  a  major  concern  must  be  directed  to  the  equilibrium  population 
established  for  the  metastable  ( 10001 )  and  (10002)  level,  which  locates 
approximately  1300  cm  1  above  the  vibrational  ground  level.  Since  there  are 
no  optical  transitions  allowed  between  those  and  the  ground  level,  a  primary 
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candidate  for  pumping  these  levels  is  the  radiationless  collisional  process 
which  occurs  only  in  the  lower  atmosphere.  The  SPIRE  data  indicate  that  the 
emission  by  the  (10001-01101 )  and  (10002-01101)  transitions  is  definitely 
observable  at  low  tangent  height  in  both  the  day  and  night  conditions .  At 
high  tangent  height  both  of  the  emissions  become  less  distinctive.  It  may 
be  concluded  that  the  radiative  pumping  mechanism  for  these  metastable 
states  is  rather  minor.  The  observation  coincides  with  the  results  obtained 
by  the  Degges '  radiation  balance  computation  program. 

The  limb-view  radiative  levels  of  the  observable  CO^  bands  for  high- 
tangent  height  are  calculated  with  a  reasonable  accuracy  by  the  model 
developed  above.  There  are  two  clearly  observable  C0^  emissions  in  the 
SPIRE  data,  one  for  the  ( 01101-00001 )  band  and  another  for  the  ( 00011-00001 ) . 
As  mentioned  above,  the  emissions  by  the  ( J. 00 01 -Oil 01 )  and  the  (10002-01101) 
transitions  are  observable  in  the  low-tangent-height  data.  Both  the 
( 10011-00001 )  and  the  (10012-00001)  band  will  be  discussed  together  with 
the  scattering  continuum  observable  in  the  3000  cm  1  ^  10,000  cm  1  region. 

The  radiance  levels  shown  in  Figures  2  and  3  are  extrapolated  from  the 

SPIRE  data  by  integrating  them  over  the  entire  band.  At  a  tangent  height  of 
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100  km,  the  ( 01101-00001 )  emission  is  3  x  10  (watt/cm  sterad).  Using  the 


expression  given  in  Eq.  (22): 


hvSO 
o  i 
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our  calculation  yields  2.0  x  10  watt/cm  sterad.  In  deriving  the  value. 


we  take 


hv  =  1.32  x  10-20  Joule 

-8  -1  2  -1 
S  =  8.3  x  10  mol  /cm  cm 

7  9 

N  =  5  x  10  x  10  for  the  100  km  tangent  height,  and 


•v  *  \ 


N  =  3  x  x  Ao  x  n  =  3  x  10^  x  .001  x  1*0, 

P 

the  Doppler  width  Ao  calculated  by 

/qvt  _t  _•> 

Ao  =  o  — ~  -  0.001  cm  for  o  =  670  cm  , 

o  d  o 

me 

and  U0  rotational  lines  as  those  effectively  emitted. 

The  SPIRE  data  indicate  that  the  ( 00011-00001 )  emission  at  100  km  tangent 

“T  2  8 

height  is  3  x  10  (watt/cm  sterad)  for  tne  daytime  level  and  3  x  10 

2 

(watt/cm  sterad)  for  the  nighttime  level.  Our  simple-minded  calculation 

_7  2  _Q 

shows  that  B  =  2.1  x  10  (watt /'cm  sterad)  for  the  day  data  and  2.1  x  10 

p 

(watt/cm  sterad)  for  the  night  data.  In  deriving  these  figures,  we  take 

_?n 

hv  =  ^  x  10  Joule 

S  =  9.6  x  10  77  mol  Vcm^  cm  ^ 

N  =  5  x  107  x  109 
o 

=  3  x  1017  x  .1  for  the  day  and  3  x  1011  x  .1  for  the  night. 

The  observed  and  the  calculated  data  for  the  day  level  agree  reasonably  well. 

The  observed  nighttime  radiance  level  is  about  one  order  of  magnitude  higher 

than  our  calculation.  In  Figure  3,  a  dotted  curve  shows  the  nighttime  radiance 
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level  calculated  by  our  formulation,  Eq.  (22).  Kumer  et  al  analyzed  the 
nighttime  zenith  radiance  data  of  the  same  band  ( 00011-00001 ) ,  which  is  shown 
in  the  same  Figure.  Both  nighttime  calculations  are  remarkably  similar. 

The  SPIRE  data  differ  from  the  data  used  in  the  Kumer1 s  analysis  in  that  they 
are  for  the  limb  view.  The  total  CO2  molecules  involved  in  both  data  differ 
by  an  order  of  a  magnitude.  Thus,  the  radiance  level  differs  by  the  same 
magnitude.  The  deviation  from  the  exponential  decay  line  in  the  radiance 
level  occurs  approximately  at  the  same  altitude  for  both  data.  Kumer  et  al 
interpreted  this  anomaly  in  terms  of  the  extra  pumping  produced  by  the 
vibrationally  excited  molecules.  The  SPIRE  data  indicates  that  the 
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difference  of  the  nighttime  radiance  level  from  the  earthshine  figure  is 

more  pronounced  in  the  high  altitude  region  where  the  interaction  to  the 

neighbor  molecules  should  decrease.  The  anomaly  in  the  nighttime  radiance 

level  for  this  band  remains  mysterious.  We  studied  a  possibility  of 

3  3 

observing  the  emission  by  the  W  Au  -  B  irg  transition  of  which  could  be 

more  intense  at  the  altitude  above  100  km  because  of  a  higher  electron 

density  there.  A  radiance  level  for  the  transition  estimated  on  a  basis  of 

the  known  data  fails  to  support  the  possibility. 

Since  at  a  low  tangent  height  the  band  is  completely  saturated  and  the 

line  structure  is  lost,  the  photon  density  available  for  the  absorption  is 

no  longer  adjusted  by  the  effective  number  of  the  rotational  lines  excited 

in  the  band  times  their  width.  It  is  given  by  the  width  of  the  entire  band 
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times  the  photon  density:  200  x  3  x  10  photons/cm  /sec.  Taking  the  photon 

—20 

energy  hv  =  h.6  x  10  Joule,  we  can  estimate  that  the  saturation  level  is 
given  by 

N  hv  c  ? 

— ^ —  =  0.9  x  10  (watt/cm  sterad). 


.^s  ..  .**  W 


Scattering  Radiance  Observable  Between  6000  cm  1  and  7500  cm-'1'  and  the 


COg-H^O  Band  at  3600  cm-1 

The  radiance  level  in  a  spectral  region  above  U000  cm-1  is  strongly 
influenced  by  the  scattering  process,  in  particular  the  level  observable 
at  a  low  tangent  height.  For  the  Rayleigh  process  the  scattering 
differential  cross-section  o(<f>)  is  given  by  Eq.  (12): 

2  2 

°m(<)>)  =  9N2°  (^2Ti)  sin2  *  3  sin2  <f>, 

where  0  is  the  optical  frequency  in  cm  \  both  n  the  refractive  index  of 

the  atmosphere  and  N  the  number  density  of  the  air  molecule  are  those  values 

at  the  S.T.P.  condition,  and  <j>  is  the  scattering  angle.  The  factor  is 

a  slow  varying  function  of  a.  The  scattering  cross-section  is  thus 

predominantly  controlled  by  0  and  <|>.  Figure  U  shows  the  front  factor 

4iro2(n-l/N)2  as  a  function  of  0.  The  radiance  level  obtained  by  the  SPIRE 

experiment  is  compared  with  the  computed  value  for  the  Rayleigh  scattering 

in  Figure  5.  The  observed  values  are  replotted  as  a  function  of  0  and  of  a 
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tangent  height  in  a  unit  of  watt/cm  strad  cm  .  Irregular-  structures 

which  can  be  considered  as  the  molecular  absorption  or  emission  are  removed 

by  smoothing  out.  In  the  computed  values,  N  is  assumed  by  the  solar  value, 
13  3  -1 

3  x  10  photons/cm  /S/cm  .  The  observed  radiative  level  supports  the 
computed  molecular  scattering  level  even  though  a  minor  difference  exists, 
dependent  on  the  tangent  height  and  the  spectral  frequency. 

The  comparison  shown  in  Figure  5  leads  us  to  conclude  two  characteristic 
features  of  the  radiance  level  in  the  3000  'v  7500  cm-1  range:  (l)  the 
scattering  radiance  is  controlled  by  the  single-scattering  process;  and 
(2)  its  main  mechanism  is  the  Rayleigh  process.  The  aerosol  particles 
contribute  to  the  over-all  scattering  process,  not  significantly  even  at  a 
low  tangent  height  of  15  km. 
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A  problem  which  we  will  discuss  below  plays  a  central  importance  in 
calculation  of  the  infrared  radiance  level  in  the  upper  atmosphere.  It  is 
a  combined  effect  of  the  molecular  absorption-emission  process  and  the 
scattering  process.  It  was  treated  rather  lightly,  if  not  completely  side¬ 
stepped.  The  computation  modeling  of  the  atmospheric  radiance  developed 
to  cover  a  region  where  only  one  of  these  processes  is  important  in  the 
radiative  process.  There  were  no  models  developed  to  treat  the  situation 
where  the  processes  are  competitive.  The  data  obtained  by  the  SPIRE 
experiment  is  in  fact  the  first  to  observe  the  radiance  level  generated  by 
these  processes  in  competition.  The  data  revealed  an  extremely  remarkable 
feature  for  the  radiance  level  in  the  region  of  the  H^O-CO^  2.7 y  band 
(3600  'v.  38OO  cm  ^ ) .  The  radiative  level  observed  in  the  15  km  tangent 
height  data  shows  an  absorptive  feature  for  the  band,  while  that  in  other 

data  shows  an  emissive  feature.  The  radiance  level  remains  rather 

—7  2  -11  2  -1 

stationary  at  10  (watt/cm  strad  y)  or  7.5  x  10  (watt/cm  strad  cm  ) 

over  a  wide  range  of  the  tangent  height  from  15  km  to  Ul  km. 

By  means  of  the  molecular  absorption  process  the  molecule  is  excited 
to  a  higher  vibrational  state  in  interacting  with  the  radiation  field.  The 
energy  absorbed  by  the  molecule  in  the  excitation  to  a  higher  vibrational 
state  is  equal  to  that  supplied  by  the  photon  field.  During  the  de-excitation 
stage,  the  same  amount  of  energy  is  released.  Dependent  on  the  environmental 
surround  to  the  molecule,  the  energy  release  takes  place  to  the  photon  field 
or  to  others.  If  the  interaction  with  the  photon  field  dominates  the  de¬ 
excitation  process,  there  is  no  net  energy  loss  in  the  photon  ensemble. 

The  absorbed  energy  from  the  photon  field  is  released  back  to  the  photon 
field  by  the  de-excitation  process.  When  the  collisional  process  dominates 
the  de-excitation  process,  the  energy  absorbed  from  the  photon  ensemble  is 


not  released  back.  The  photon  density  in  the  ensemble  decreases  until  the 


photon  ensemble  reaches  a  thermal  equilibrium  with  the  molecular  ensemble. 
The  photon  density  which  is  established  under  a  thermal -equilibrium  is  given 
by  the  expression  which  was  used  to  compute  the  values  for  Figure  1: 


2tto2c 

N  =  — - -  . 

p  hoc 

kT  , 
e  -1 

The  radiative  level  observable  from  such  an  ensemble  of  the  photons  and 
the  molecules  is  given  by 


_  hvN 

1  =  _ £  • 

IT 

The  stationary  emission  observed  for  the  2.7p  H^O-CO^  band  corresponds 
to  the  situation  described  above.  What  we  can  observe  is  the  photons 
generated  by  the  H^O-CO^  band  which  are  in  a  thermal  equilibrium  at  the 

7 

tropopause  temperature.  Using  the  blackbody  photon  density  N  .92  x  10 

P 

—1  —20 
at  3600  cm  *>nd  225°K,  and  the  photon  energy  hv  7.1  x  10  Joule  at 

—1  —11  2 
3600  cm  ,  the  saturation  level  I  is  calculated  as  2.1  x  10  watt/cm 

sterad  cm  With  assuming  T  =  260°K,  the  computation  yields  5  x  10 
2  -1 

watt/cm  sterad  cm  .  This  value  is  independently  determined  from  the 
internal  excitation-de-excitation  mechanism  of  the  molecular  system. 

Whenever  the  molecular  system  establishes  a  thermal  equilibrium  with  the 
environment  under  the  saturation  absorption  condition,  the  photon  field  must 
be  in  an  equilibrium  with  it.  The  stationary  photon  flux  observed  for  the 
°.7y  H^O-CO^  band  represents  this  situation,  as  was  the  case  for  the  CC>2 
( 00011-00001 )  transition  at  lower  altitude.  If  the  saturation  photon  flux 
is  larger  than  the  scattering  photon  flux,  the  radiance  level  of  the  band 
in  question  is  emissive.  If  it  is  smaller,  the  radiance  level  is  absorptive. 
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In  both  cases,  the  radiation  level  remains  constant. 

The  Radiance  Level  of  a  Molecular  Band 

In  the  last  two  chapters,  we  discussed  the  saturated  radiance  level  of 

the  band  of  CO^  and  of  the  2.7u  H^O-CO^  band.  We  found  that  the  photon 

flux  under  the  saturation  conditions  is  given  by  the  blackbody  radiation 

formula;  i.e.,  that  the  molecular  system  and  the  photon  ensemble  come  to 

establish  a  thermal  equilibrium.  For  a  given  molecular  band,  there  exists 

a  critical  altitude  that  the  photon  density  is  completely  specified  by  the 

blackbody  radiation  formula  independent  from  the  vibrational-rotational 

transition.  In  the  atmosphere  above  the  critical  altitude,  the  photon 

density  is  specified  by  the  molecular  transition  equation,  Eq.  (U): 

N  =  It  H  N  S. 
t  os 

A  sensor  placedvery  high  in  the  atmosphere  collects  the  photons  generated 
by  the  molecular  transitions  which  occur  along  a  line-of-sight .  With  the 
approximation  of  a  thin  optical  density  air,  the  total  molecular  transitions 
which  contribute  to  the  observation  are  given  by  integrating  the  transitions 
along  the  path: 

N  =  /  N  dfc  .  (23) 

Under  the  approximate  expression  derived  for  an  observation  at  a  tangent 
height  x,  we  calculate  the  total  density  per  a  solid  angle  per  a  unit  cross- 
section  by 

1*N  S 

N(x)  =  — 2-  N  (x)  2  /2H( x+p )  .  (24) 

IT  O 

In  deriving  the  equation  above,  an  assumption  is  made  that  the  photon  density 

N  responsible  for  exciting  the  molecular  transition  is  constant  along  the 
P 

path.  In  a  spectral  region  below  2000  cm  ^  where  the  telluric  photons  at 


the  tropopause  temperature  control  the  excitation,  the  assumption  is  obviously 


well  justified.  In  a  higher  cm  ^  region  where  the  solar  photons  are 
dominant,  a  blackbody  temperature  corresponding  to  the  solar  photon  density 
is  much  higher  than  the  tropopause  temperature  as  seen  in  Figure  6.  The 
radiance  level  of  the  molecular  emission  reaches  a  maximum  plateau  at  the 
saturation  photon  density  specified  by  the  temperature  of  tropopause.  Thus 
the  molecular  transitions,  even  when  they  are  integrated  along  a  long  path, 
are  much  less  than  the  solar  photon  density  flux  which  has  a  temperature 
much  higher  than  the  saturation  photon  temperature.  The  loss  caused  by  the 
molecular  transition  along  the  path  is  therefore  much  smaller  than  the  total 
solar  photon  density,  and  a  constant  photon  density  is  maintained,  as  assumed 
above.  This  justification  obviously  breaks  down  as  the  radiance  level 
approaches  the  saturation  level. 

In  a  non- saturating  region,  we  find  that  the  observable  photons  given  by 

Eq.  (22)  are  controlled  by  the  factor  N(x  ),  the  molecular  density  at  the 

o 

tangent  height  X.  In  the  first-order  approximation,  the  radiance  level  at  a 
particular  tangent  height  gives  the  density  of  the  observing  molecule  at  that 
altitude.  For  a  stable  molecule,  its  density  decreases  exponentially  as  a 
function  of  the  altitude.  We  can  model  the  radiance  of  a  molecular  band,  as 
snown  in  Figure  T,  by  an  exponential  function  with  a  plateau  at  the  critical 
altitude.  The  model  certainly  agrees  with  the  radiance  data  collected  by  the 
SPIRE  for  the  v,  band  of  CO^.*  The  computed  results  by  the  Degges  program  as 
well  as  by  the  LOWTRAN  for  various  bands  support  this  model. 


*The  atmospheric  sounding  technique  currently  experimented  uses  this  model. 

It  is  to  observe  the  saturation  photon  density  at  the  critical  altitude,  which 
varies  from  a  spectral  region  to  another  as  a  function  of  the  absorption 
strength.  The  observed  radiance  level  gives  a  local  equilibrium  temperature 
at  the  altitude.  An  atmospheric  sounding  program  proposed  to  NESS  will  use 
the  (01101-00001 )  and  ( 00011-00001 )  transitions  of  CO^.^ 
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Figure  7 
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In  the  upper  atmosphere,  the  radiance  level  is  thus  dictated  directly 
by  the  molecular  density.  A  deviation  of  the  radiance  level  from  the 
exponential  decay  is  expected  if  the  molecular  concentration  falls  off 
rather  differently  from  the  exponential  model.  The  SPIRE  data  on  the  infra¬ 
red  radiance  level  provide  a  base  to  determine  the  concentration  of  the 
observed  atmospheric  species  as  a  function  of  the  altitude.  From  our  point 
of  view,  the  radiance  level  data  of  0^,  H^O  and  NO  and  other  minor  species 
provided  the  first  measured  data  which  revealed  the  distribution  of  these 
species  at  high  altitude. 

The  radiance  level  data  of  0^  observed  at  a  tangent  height  of  80  km 
show  a  rather  dramatic  change  from  the  day  to  the  night  side,  while  the 
density  of  the  excitation  photons  remains  unchanged  between  these  two 
periods.  A  mechanism  which  causes  the  diurnal  change  must  be  sought  in 
some  photo-chemical  effect  triggered  by  the  ultraviolet  photons. 

Modification  Applied  to  the  BGND  Program 

In  the  course  of  our  study,  we  made  an  extensive  use  of  four  computation 

11,12 

programs:  the  BGND  written  by  Degges  of  Visidyne,  Inc.;  the  LOWTRAN  by  Selby 

l1*  5,6,7 

of  then  AFGL;  the  FLASH  by  BlKttner  of  Radiation  Research  Associates,  Inc.; 

15,16 

and  the  FASCODE  by  Clough  et  al  of  AFGL.  We  applied  some  modifications  to 
these  programs  for  improvement  of  the  execution  efficiency.  The  operating 
system  of  our  central  site  computing  facility  is  designed  to  accommodate  a 
large  number  of  the  remote  terminals.  Because  of  the  operating  system,  we 
found  that  the  program  executes  more  efficiently  if  its  working  space  required 
is  reduced  in  size  from  the  original.  We  divided  these  programs  into  small 
sections,  achieving  their  executions  through  our  remote  terminal  possible 
during  the  regular  hour.  In  addition,  inefficient  10  operations  found  in 
these  programs  are  streamlined  to  meet  our  specific  needs  and  to  reduce  a 
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computational  waste.  We  have  consequently  achieved  a  substantial  improvement 
in  running  these  programs.  Our  strategy  for  using  the  FLASH  program  was 
described  in  our  Report  No.  1.  Overall  modifications,  except  for  the  BGND, 
did  not  extend  to  rewrite  their  computation  algorithms.  We  feel  that  they 
have  no  specific  merit  for  their  coverage  in  this  report.  The  case  for  the 
BGND  is  different.  Our  modification  is  listed  in  the  Appendix.  A  general 
logic  flow  of  the  program  execution  is  illustrated  in  Figure  8. 

Other  Efforts 

(1)  For  the  balloon-borne  experiment,  we  evaluated  the  test  data  on  the 
interferogram  PCM  telemetry  scheme.  We  found  some  mysterious  bit 
droppings  in  the  PCM  interferogram  data  string.  A  substantial  effort 
was  poured  to  resolve  this  mystery  in  a  joint  effort  with  the  Idealab 
and  the  AFGL.  No  conclusions  were  made  on  this  problem.  Our  effort 
was  switched  after  a  new  contract  for  the  balloon-borne  data  analysis 
became  effective. 

(2)  A  very  efficient  CRT  display  program  for  the  CYBER  system  was  written 
by  R.  Lauzzana.  Its  listing  is  given  in  Appendix  B. 

(3)  In  anticipating  a  large  amount  of  the  CYBER  generated  plots,  we 
implemented  a  micro-computer  controlled  linkage  between  our  plotting 
machine  and  the  CYBER  MODEM  scheme .  The  program  implemented  on  ROM  of 
the  8080  system  is  listed  in  Appendix  C. 

Conclusion 

The  SPIRE  experiment  is  in  a  sense  a  historical  event  in  the  progress  of 
atmospheric  physics.  Even  though  the  data  were  collected  using  a  simple 
spectrometer  radiometer,  their  merit  is  extremely  valuable  in  developing  our  under¬ 
standing  of  the  atmospheric  radiative  process.  Our  simple-minded  model 
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Figure  8 


achieved  a  good  success  in  understanding  the  radiance  level  of  the  major 


atmospheric  CO^  hands.  We  believe  that  our  constructed  model  provides  a 
good  insight  into  the  atmospheric  radiative  process.  Ours.,  is  no  different 
from  the  aforementioned  computation  programs ,  since  all  of  these  are 
constructed  on  the  same  theoretical  develojment .  Ours  attempted  to  model 
the  radiative  process  in  the  simplest  form  possible  required  for  providing 
a  good  estimate  of  its  level. 

There  will  be  no  question  that  the  data  will  be  collected  with  an 
improved  spectral  resolution,  as  time  will  progress.  Then  we  will  come  to 
understand  more  details  of  the  radiative  process  in  the  atmosphere.  In  so 
doing,  we  will  be  more  appreciative  of  those  aforementioned  computation 
programs.  At  this  moment,  our  simple-minded  model  functioned  reasonably 
well  for  analysis  of  the  SPIRE  data,  revealing  a  basic  character  of  the 
infrared  radiance  level  in  the  upper  atmosphere,  as  well  as  some  abnormal 
feature  which  must  be  studied  more  carefully  in  the  future. 

The  following  list  summarizes  our  model  of  the  radiative  process 
in  the  upper  atmosphere: 

(1)  In  the  upper  atmosphere,  the  molecular  absorption  transitions  which 
trigger  the  emission  are  well  approximated  by  assuming  that  they 
originate  from  the  ground  state. 

(2)  The  photons  available  for  exciting  the  molecular  transitions  are 
either  those  of  the  telluric  origin  or  of  the  solar  origin.  The 
cross-over  takes  place  in  the  vicinity  of  2000  cm~\ 

(3)  Below  the  critical  altitude,  the  molecular  bands  establish  a  local 
thermal  equilibrium  with  the  photon  field.  The  radiative  level  given 
by  the  blackbody  radiation  formula  is  maintained  over  a  large  range 
of  altitude  if  the  absorption  is  very  strong.  It  is  independent 
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of  the  observation  condition  and  of  the  scattering  radiance  level. 

(U)  Well  above  the  critical  altitude,  the  radiative  level  of  the  molecular 
band  is  proportional  to  the  molecular  concentration  at  that  altitude 
in  the  first-order  approximation. 

(5)  The  scattering  radiance  level  in  the  near  infrared  region  is  primarily 
determined  by  the  Rayleigh  process  with  a  single  scattering  center. 


References 


1.  A.T.  Stair,  Jr.,  R.  Nadile,  G.  Frodsham,  D.  Baker,  and  W.  Grieder, 

Atmospheric  Emission  Spectra  (l.it3-l6pm)  of  the  Earth's  Limb  (SPIRE), 
Topical  Meeting  on  Atmospheric  Spectroscopy,  OSA  Proc.  (1978). 

2.  P.  Hansen,  R.  Lauzzana,  H.  Sakai,  and  J.  Strong,  AFGL-TR-78-0013  (1978). 

3.  R.M.  Goody,  Atmospheric  Radiation  I,  Oxford  (1964). 

4.  E.J.  McCartney,  Optics  of  the  Atmosphere,  Wiley  &  Sons  (1976). 

5.  Blattner  et  al,  Appl.  Opt.  13,  534  (1974). 

6.  Blattner  et  al,  AFCRL  Report  TR-75-0317. 

7.  Blattner,  RRA  Research  Note  RRA-N7512. 

8.  Pauling  and  Wilson,  Introduction  to  Quantum  Mechanics,  McGraw-Hill, 

New  York  (1935) . 

9.  S.S.  Penner,  Quantitative  Molecular  Spectroscopy  and  Emissivity, 

Addison-Wesley  (1959). 

10.  The  transition  notation  for  CO^  is  conformed  with  the  AFGL  Atmospheric 

Absorption  Line  Listing. 

11.  T.  Degges,  AFCRL  Report  TR-7U-O606. 

12.  T.  Degges  and  H.  Smith,  AFGL  Report  TR-77-0271. 

13.  A.  Kaplan  et  al,  Appl.  Opt.  16,  322  (1977). 

lit.  J.  Selby  and  R.  McClatchey,  AFCRL  Report  TR-75-0255. 

15.  S.A.  Clough  et  al,  AFGL  Report  TR-77-0164. 

16.  S.A.  Clough  et  al,  AFGL  Report  TR-78-OO81. 

17.  J.  Kumer  et  al,  J.  Geophys.  Res.  83,  4743  (1978). 


32 


sendix  A 


PROGRAM  EGREACCOLTPU'T.TAPEl.T  APE2.TA  FE6) 

THIS  IS  THE  FIRST  PART  CF  CEGG*S  fGNC  PROGRAM.  IT  JUST  REAOS  CEGG 

cats  Inc  gei'S  oata  react  for  frogram  qanrad. 

CAST  CHANGE  AUG  10,  197* 

C IHENSIC  N  XIC(4,25),XICT (4) ,  A  IT  C  20)  ,  STORE  (  70  >  ,  IAL  T  (2  5  > 

DIMENSION  CAT  ATT  <2$) ,CC  EX  T  (EE)  .  CCEX9  (25) , NLEVEL ( 2 5 > , NP ANDt 25> 
DIMENSION  LSC 125, Hi  ,LCeC25,2!),FCM(25,25)  ,8 VC 25, 25)  ,AMASS(25) 
CIHENSICN  2ANGC7I  ,STAU<?) ,SFAC  C81.RCN  C8) , RATI CCA) 

C| PENSION  LECC25 , 25)  ,L  El  (  25 , 2  5  )  , STR  ( 25 ,25)  ,SFL IX C 25,  25 ) , TFLUXT 25 , 

CIHENSICN  CC2(18,2),WATIA<1?,2) 

CIMENSION  CHEPLCiei) 

REAL  LAM  C A 

COMMON  /VeRCCM/  I  ,CC2  ,  K ATv A , CCEXT, C CEXR 
COMMON  /TEQCCM /  NLTE 
CCMMCN/K  CEM/KCHEM 
CCMMCN/I NCUT/NR ,NF,NN 
1^MM^N<C^NC  C201, 8),T<201),LAMCAC?5>  , 

CATA(ZANC(I)»I=1»7)/C.C,<0.C,4C.G, 60.0,70.0,  80.0,90.0/ 

GATACCAT  ATYCI),I  =  1,25) /ENMr THANE  ,8HCAPBCN  0,  f NHATER  V A , 8HN ITR  TC  C 
1 ,  EHNlTRO  US  ,?H07CNE  PR.fFCZCNE  NC, 18*81-  / 

KCHEN  =  1 
NR  =  1 
NP  =  7 
NT  X  =  6 
NW  =  6 

10  FORMAT  C 12) 

READ  C  NR.  1C)  IGAS 

20  FCPMATC1H.22HNUKEER  CF  GAS  FFCFILES,  13) 

WRITE  CNN  .20 )  IGAS 
IGAS=IGAS41 
GO  17  C  J  =  1 »IG AS 
30  FORMAT  C 4 AE . 1 3. F7  ,  C ) 

READCNP, 20)  CXIDCK,j)  ,  K  =  l,4),  TCHK.ALTM 
IF (J.EQ.ICHK)  GO  TC  50 

40  FORMAT  C  7HC  CAPO,  A6  ,  2  OH  CUT  CF  CRCER  FOR  GAS,  13) 

WRITE  CNN  ,40MXIOCN,4)  ,K  =  1,4)  ,  J 
50  CC  105  K=l,€ 

K1=(K-1) *541 
K2=K1 44 

60  FORMAT  C5CF5.0,F8,1,2X) ,12,13) 

READCNR, E0> CALT  CL),STCRE CL)  ,L=K1,K2)  ,  ICHK,  KCFK 
IFCICNK. EC.J)  GO  TC  80 

70  ££PMAT(3£HCATA  CARC  CUT  CF  CPCER  FCR  GAS, 13,  3X,  BHCARC  NO.,  12) 

KF1TE  CNN, 70)  J,IFIVE 
80  IFCKCHK. EC.K)  GO  TC  50 
WRITE  CNW  ,70)  j, K 

iFcKtcLKI^ALTM)  GC  TC  99 
KSET -L 
GC  TO  120 
99  CONTINUE 
100  CONTINUE 
IDE  CONTINUE 

110  FCRMATT4 <N*AX  ALT,  NOT  FCUNC  ON  DATA  CARDS  -  GAS  NC.  ,  12) 

WRITE  (NW  ,110)  J 

Ho  £$0&E°(lV  =iCcG*ST  CFE  Cl)  ) 

IFCALTM.GT.19C.)  GC  TC  140 
JALTC j)  =  AlTH-59. 

GC  TC  150 

140  IALTC J)=Sl.4f ALTM-15C.) /5. 

150  CALL  INTERFCALT, STORE, KSET, l,v) 

.  tHsm.i. 

160  CCNC ( l,J  )  =  EXF CCCNC CL,J ) ) 

170  CONTINUE 

CC  175  J  =1,201 
17!  TCTOENCJ  )=CCNCC J,  IGAS) 

R  E  AC ( NR, 20)  CXIOTCK),K  =  1,4),  ICHK,  KCHK 
IFCICNK. EC. 1)  GC  TC  ISO 

1*0  FCRMATC  7HIC  CARC,  4A6,2€HCLT  OF  ORDER  FCR  TENF  OATA) 

190  U yUW'A'V  »«^K),K=l,4) 

K 1= ( K-l)  *541 
K2=K1  44 

ffm^:i5,.i,raLUsJ5r(U’L'K1’K?l’  ICHK’ KCH< 

200  FORMAT  C  3 1HTEMP  DATA  CAFC  CUT  CF  CRGPR  NO.-th  33 


GAS  NC.  ,  12) 


y  v’v 


no  HHHRiHSI  g8  tc  25C 

2  2  0  WRITE  (NW,2fi«)  K 
230  DC  240  L=K 1 ,K2 

JFIAL71L  J.LT.ALTH)  GC  1C  240 

gPIc126  C 
240  CCNTINUE 

ft!  FCRMllVf 7HHAX  ALT.  KCT  FCUNC  CN  TEMP  OATA  CARCSI 
WRITE  INK  ,250) 

260  CALL  INTERF<ALT.STCRE,XSE7,2,1) 

261  FORMAT  ( IX ,15, I F «E 12. 4  ) 

WRITE  (NW  ,261)  C  J  ,  (CCNC  <  J  ,  I )  ,  1=  1 ,  IG  AS )  ,  T(J),J  *  1,201) 

IGAS=IGAS-1 
DC  331  11*1, IGAS 

CC  311  J=l,25 

IFfTYFE. EC.CATATY  (J) )  GC  TC  312 

311  CCNTINUE 

370  Fc£mAT(32fTHE  GAS  IS  NCT  IISTEC  IN  THE  PROGRAM) 

STOP  1 

312  NSF=J 
NLEVEKN  SF )  =  NLEV 
NE  ANC  (NS  P) =  NPAN 
AN  ASS (NS  F)  =  AN  AS 
CCEXT <NSFI=CEXT 

?^NiF^icl25^AD(KR,  320  )  (TC?  (1 ,  1 ) ,  C  C2  (1 , 2 ) ,  I  =  1,  1  8) 

320  FCRMA1  (6CF5.2,F8  ,2)) 

IF (NS F. EG. 31  REA 0<  NR,  32  0  )  (W A  TV A  (1,1)  ,WATVA (I , 2 ) ,1 =1, 18 ) 

IIFVL  =  0 

CC  313  T  »=  1  ,NLEV 

READUR^  33  0)LSC(  Iu»N  SF  )  ,LCC  (I  J  ,NSP)  ,  RC  M IJ  ,N?F),9Vd  J,NSP>  ,  JGA  J  * 

230 1 ( 2IE,2E12.4,36X,2I3) 

IF (NSF.EC  .jtASI  GC  TC  314 

340  FCRMAT(64HTHE  COCE  CN  TEE  OATA  CCES  NOT  AGREE  WITH  THE  FIRST  CARO 
1WHICH  IS, 13) 

STOP  2 

314  IFdLEVL  .EC.JLEVEL)  GC  TC  313 

350  FCRMAT(42HTH£REN!S  A  CARC  CUT  CF  SEQLENCE  IN  GAS  NC.  ,13) 

STCP  3 

313  CCNTINUE 
IlEVL  =  0 

CC  316  IJ=1,NEAN 

RE  AO (NR^  360 )LeC( I J,NSP ) ,L  EU (IJ ,NSP>  ,  STR  (I  J  ,N  SF ) ,  SFLU  X  (I  J  ,NSP) , 

367  FORMAT  ( 2T 6 , 1PE 1 2 . 4 ,21 3 , IX , 4 A6) 

,c  5" 

STCP  4 

222  IFdLEVL  .EC.JLEVEL)  GC  TC  316 
WRITE  (N  W  ,350)  NS F 
STCP  5 

316  CCNTINUE 
331  CCNTINUE 

371  FCRMAT  (1F6E17.4) 

READ (NR, 3781  CHEML  „ 

30  REAO (NR, *)  IGAS, NLTE 
IF (IGAS. EO. C)  STOF 

KJHISSjftH  T  O  T  C  E  N 

WRITE (NP ,376)  ( CCNC (NZ, IGAS) ,N7= 1, 20 1) 


ecu  ILCOl  JL  winjr  l  u  LU  uuiiuri  t  ~ 


WRITE  (NP  ,37  3  )  ( CCNC (NZ, IGAS1 ,NZ= 1, 20 1) 

WRITE  (NP  ,360  IGAS, NLTE 
NLE  V  =  NLEVELdGAS) 

WRITE  (NP^I’O  <(XlC(K,IGAS),KSl,4),NLEV,NPNC,AHASSdGAS) 

2  ,CCFXT(  IGAS), COEXEf  IGAS) 

FCRMAT(4  A€,6X,T2, 3X,I2,3X,F6. 2,4 X, 2<E12.4,3X> ) 

°WR?TF  (NF^eJ)  Ktl)fdJ,IGAS),LCC(IJ,  IG  AS)  ,RCM  (IJ  ,  IGAS )  ,EV  d  J,IG  AS ) 
CCNTINUE 


C  382  Ij  =  1,  N  BFC 
382  CONTINUE 


* IGAS ) ♦ STR 1 IJ* IGA  S) 


i  J  A  IT 


IF  ( 

GEITE  16 »  378 )  CHEMl 
GC  TO  83  0 

ENC 

StBROlTIFE  INTER P  t Al T , E TORE *K S FT ,MOD E» I) 

y^cFScNcHSl  Jiffl^SHSiAh-CA  (25)  t 
1TCTDEM2  Cl) 

J  =  1 
A=60 . 

8  CC  120  K=1,K$ET 
10  IF  (A-ALT  (lOieC,  20  ,20 
20  GC  TC  (3  C«4C»  ,MOCE 
30  >  =  STCRECK» 

GO  TCJ) =STCfi£(K> 

50  IF  (K.LT. KSETI  GO  TC  60 
RETURF 

€0  IF1J.CE.E1)  GC  TC  70 
A  =  A  ♦  i  . 

GC  TO  UC 
70  A=A*5. 

GC  TC  lit 

80  FF=(A-ALT1K-11)  /  1  AL  T IK  )  -  ALT  <  K  -  1 » ) 

GC  TC  <9C«100),  KCCE 

90  CCNCC  J,I  )  =  £TCRE  (  K-l  I  4FR  *  <STC»E  ( Kt  -5T  CRE  1  K- 1 )  ) 

91  J=J*1 

IF(J.GT.«1)  GC  TO  E2 
A  =  A«-1 . 

92  KM.19 

GC  TC  10 

100  7  <j>=STOCE<K-1)»FFMSTCFE<K|-STCFECK-1>> 

GC  TO  91 

110  J=J4l 

120  CCNT IFUE 
RETURF 
EF  C 
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FRCGRAK  FANRAC(TAFE2  =  1C2E,TAFF4=1Q'’B,T/PE3,OUTFUT  =  102E 
2  ,TA  F  E6=  C  L  T  FL  T) 

THIS  IS  A  ^COIFIEC  VERSION  CF  CEGG'S  SUPRCL'Ttnf  3ANRAC 
LAST  CHANGE  AL'G  3C,  1R7? 

G1NENSICN  FCKER  <40  ,  KL  (AC)  ,XIC(4> 

CIMENSICN  TFLlMf  (  AC)  ,?FLUxUo  I 
CINENSICN  GF  AC  C  4  0  > 

CINFNSICN  LSC  (4<J)  ,LCC(4C)  ,L°U  <4G),LeC<40) 

CIVENSIC  N  VAVE<40),STR<40t,ev<40>,QNCLAM(4n>,FACC(40) 

Cf^FSICN  GNCSTF(2C1),TEMF  r?01  J.CENTCT  (201 >, vie(40,3  > 

CIMEN?ICN  ENCSA0(2tI  ,4C>  .BNCTAlUcI,  40 
CCMMCN  VIEFCPt 20 1 ,2!) 

CCNMCN  /  CCL  /  TVSI0(2C1)  ,VTSIG  (201) .VVLSIG (201)  .VVCSIC  (201) 

2  .VVNSIG  <2C  1)  ."VRSIG  <2  C II 

CCMMCN  /CFC8TA/  V  ALT.ME  VEL  .NBANC,  7FAC,KP  (4  0)  ,EA  (40  )  ,CUF(A  C) 

2  ,GCCWN(40),TTUP(IC),TTCCNN(4C),LFSTA<40),  GOWN  (4  0  )  ,T  TM2C1 ,4  0) 

2  .C0NCC2 Cl) ,CHEML  (201) 

FCUIV&LENCE  (GKCETE(l)  ,VIPFCF(1,1)  ) 

ECLIVALENCE  <  TFL  L  X  ,  GLF  )  , ( SFLUX ,GCCWN  ) 

ECLIVALENCE  (TTW,ENC1AL)» (TTH.ENCRAC) 

ECLIVALENCE  <LSC(1), 71^(1,11)  ,  (LCCd  ),VIE(1,2)),  (LEU  <  1)  ,  VI  0  (  1,  ?)> 
PC  1 1  VALENCE  (GFAC.CCKN)  ,  (  CFNTCT  (  1  )  ,  T  TM  1 , 1 1 ) 

EXTERNAL  CC  2  ,  H2  C  ,  F  N  C  ,  C  2 
ZF  AC  =  1  •  0 

READ  IN  TEWFERATURE  ANC  NCIECLLAF  PRCFILE 

READ(2,37S)  TEMP 

WRITE  (4t  278  )  T^MP 

READ(2,37M  CENTCT 

READ(2,37fi)  CCNC 

FCFMAT(1F€E12.4) 

READ  IN  VVICH  GAS  ANC  CCNCITICN  ( NL  T  F=  3  NIGHT,  NL  TE=  4  CAV) 
REAO(2,368)  IGAS.NLTE 
WRITE (4,  26  P  )  IG AS  ,NL  IE 
FCFMAK2IE,  ic5c12  .4) 

NCTE  I  IF  C EGGS  IS  ICAS  IN  THIS  FPOGFAN 
READ  IN  LEVEL  TNFC 

RE  AC  (2, 3  1C)  (  XT  C I  *  >  ,  K=1  ,4)  , NL E VEL ,N P AN C, RM . CC F XT , C CE XE 
WRITE U, 21G)  (VIC  (K J  .K=i,4) , NLEVEL .NEAND.RN, CCEXT , GO  EVE 
FCRMflT(4A6,EX,I2,2X,l2,ZX,F€.Z,4X,2ClFEl2.4,?X)> 

CC  1  1  =  1  ,MEVEL 

REAC(2,3€8)  LSC( T  )  ,LCC  <T>  ,WAVF (I) ,?V (I) 

^I5E^ce?^bscu,’LCCa,’KflVE<I,’RV(I’ 

CC  2  1  =  1  ,NEANC 

EEAO  (2,3  E6  )  LEC  (I  )  .LEI  (I)  ,STR  (I)  ,SFLLX  (I)  ,  TFLLXI ) 

REAO(2,3ee)  NALT.JlT 

IF(IGAS. EC.4)  RE«?(2,5C)  (CVENL(T)»I=1,161) 

F  CRN  A  T  ( E>  E  1 2  .4  ) 

THIS  SECTICN  CCHFLTFS  FCPLLATICN  CF  VIPRATICNAL  LEVFLS. 

CCMPU7E  THERMAL  FCFLLMICN 
CC  1 3  C  N  7  =  1,  KALI 
SLM  =  G.C 
THP  s  TEPF(NZ) 

CC  11C  NL  =  p,  NLEVEL 
XVN  =  H  CC (LSCCNL )/10, 1C) 

AA  =  EXF  (-1,4387S*WAVE  (KL) /THF) *XYN 
VIBFCF(KZ,M)  =  AA 

i  tMr  •  '* 

SDH  =  SIN  ♦  1.0 
CNCSTE(NZ)  =  CCNC (NZ) /  S  L  N 
,  Z»NL)  3  CCNC(NZ)  *VIEFCP  (NZ,  NL  >/SUN 

!  t  CN  T INU  E 
)  CONTINUE 

FIND  LEV  Et  CCUPLFt  WITN  NITROGEN 
CC  1 3!  NL  =  2,  NLEVEL 

IF  ( ICC  (NL )  .Ed.  1)  VVVAVE  =  2’31.-  WAVE  (NL ) 

IF  (LCC(NL)  .EC.  2)  VVVAVE  =  1*56.4  -  WAVE(M) 

FAC  =  2  331  .0 

IF  (LC,C(NL)  .EQ.  2)  FAC  =  1556.4 
CCNTINUF 

IF(NLTF.EC.t)  GO  TC  17C 

FINO  COL  IISICN  EXCITATICN  ANC  CE -EXC IT  AT  ION  CCFFS 

cSfJ8JznI  hUlM 

TEMF13  =  FXP(ALOC(TEHFNZ)/3.D 

IF  (IGAS.EC.  3)  GC  TC  136 

IF  ?»  CO  TC  140 

VT5nIy7)  5  6.6<)E-lC*FXP(-e4.07/TEMF13) 

VVRSIG(NZ)  =  1 . 7l€-6*FXF  (-175  .3/TEMF1J) 


EQ. 

1) 

VVVAVE  =  2  "*31 .  - 

WAVE  (NL  ) 

EC. 

2) 

VVVAVE  =  1556.4 

-  WAVE(M) 

EQ. 

2) 

FAC  =  1556  .4 

GO 

TC 

1  7  C 

iFjfVSw 

TEMF13  =  FXP  (ALOC  (Tl 
IF  (IGAS.EC.  3)  GC 
IF  i J5A  1*NE .  2)  GO 
=  6 . 6R  £  -  1 
VVRSIG(NZ)  =  1.71E-I 


VRSIG(NZ)  =  1 . 71€-€*FXF  (-175  .3/TEMF1J) 
♦  €.07E-14*EXP(  15.27/TENP13) 


YVDSIktKZI  *  l.Ot-li  4  ‘.16t-ll*tXH  <-7b.75/TfcKPl-JI 

GO  TC  150 
CONTINUE. 

VVNSIcJkI)  =  4.  CE-1Z 

VTSICCN 2)  -  5.376*10*  EXP(-70  *  0 /TEMP  13) 

^R^GCK|I  =  VVNSIGTKZl’EXPl-l.MSaZS/TEMR  WZ)  *1594.  736) 

CCNT  INul 

VTSI G< N  2)  =  CDEXT 
).VNSIG(N7)  =  CDE  XE 

^RFJGJNZ,  =  VVNSIG(N2)*EXP(-1.43879/TENPfNZ)*VVWAVE) 

TVSIG(NZ)  =  VTSIG  (NZ) /GNCSTE (NZ )*VI EPCP (NZ, 2) 

VTSIG(NZ)  =  VTSTG<NZ)*CENTGT(NZ) 

TVS  I G (N  2 )  =  TVSTG(NZ) *OENTCT(NZ) 

TNNV  =  TEMF(NZ) 

IF  (NLTE  .IT.  51  GO  TC  154 

IF  (NZ  .GT .  41  .AND.  N2  .IT.  51)  TNNV  =  TEMP  (N7  )  4FLCAT  (NZ-41 ) 
L  *280. C 

IF  (NZ  .GE.  51)  TNNV  =  3CC0.C 
CONTINUE 

VVN^5G(Nzf^=*VVN^l?lNZ)*?EN7CT(NZ)*FACTOR/(l.C  ♦  FACTOR) 

V VRS IG ( N  2 )  =  VVRSIG(N2)*CFNTCT(N7)/ (1.0  4  FACTOR) 

CONTINUE 

COMPUTE  VIEFCP  FCF  7  EMF  =  1 8  0 
DC  166  N  2  *  1,  NAIT 
SUM  =  0.0 
7  MP  =  1 f 0.0 
CC  162  NL  =  2.  NLEVEL 


VIEFCP  FCF  T EMF 1 1 8 0 
2*1,  NAIT 


(NL) /TNP) *XYN 


SUM  =  SIM  4  AA 
162  CONTINUE 

SUM  =  SIN  4  1.0 
GNDSTE(NZ)  =  CCN((NZ)/SUM 
DC  164  NL  =  2,  NLEVEL 

VIPFCF( NZ.NL)  =  CCNC(NZ)*VIBFCP (NZ, NL) /SUN 
164  CONTINUE 
16 (  CONTINUE 
170  CONTINUE 

?8h?wen^i5;in?  AN  C 
NBC  =  LEC(NB) 

NI  =  MO  C(NEC/1  00 ,100) 

NF  =  MOC(NEC,100) 

wmi  i  ttmmvAh10' 

TNAVE  =  NAVE ( N I )  -  VAVE(NF) 

E  NOL  AM ( N  B )  =  1.0E44/TVAVE 

IF  ( TFLLX  (NE )  ,EC.  0.0  TFLLX(NB)  = 

1  5.9S7E46*TNAVE**4/(EXF(1.47C79*THAVE/TEMP( 1))  -  1. C) 

XYM  =  MCC(LSC(NF)/10,1C) 

XTMM  =  MCC  (l<jC( NI  I/1G. 1 C) 

E  A  (NE)  =  STf$(NE)*2.804/ENCLAM  <NE)**2*XYM 
1  *E XF  ( 1  .4?879*WAVE(NF1 /2B6.C)/XVMK 
FCWEF(NE)  =  1.580E5E-2C/BNClfN(NB) 


GFAC(NB)  =  3.72E“24*STF  (N'E)*ENClAM(  NB  )**? 

*EXF(1.43872*WAVf(Np)/206.C) 

CUP(NB)  =  GFAC(NEfMFU.X(NE)*  12.4864 
IF  (NLTE  .LT.  3)  GLF(NE)  =  C.C 
GCCWN(NE)  =  GFAC  (NE:*SFLUX  (NE) 

IF  (NLTE  .LT.  4)  GCOWNlNB)  =  0.0 
K8N  =  KE(NE) 


XLN  =  K  L (NE) 

E VN=BV ( 1  ) 

ENQL  =  ENCLAM(NB) 

F ACC  (N8 ) =  POWER  (Ne)*EA  (Ne) 

SNB  *  STR(NS)*EXF(1.4v?79*NAVF(NF)/296.Q) 

DC  18  C  N  2  =  1,  NAtT 

ENDTAU( NZ.NB)  =  TAUMAX(KLN,KEN»RM.8VN,TEMF(NZ),  BNCL  »SNP, 
1  VI B FOP  (N2,NF) ) 

IF  (NLTE  .GT.  1)  GO  TC  18C 
^NO^AC(NZ,NB)  =  VT€FCF(NZ,M)*FAC 

CC  175  N  E*  1  ,NEANC 

WFITF  (4, Zi *)  LBC(NE)  ,LFL  (NB)  ,£TR  (NB)  ,eNCL  AW(  Ne ) ,  F  ACC (NP) 
WRITE  (4, 36  8)  NALT.JALT 


i  e  i 


182 


190 

200 


10G0 

10*0 

1  CGI 

1 C  C  8 
210 
220 

230 

2<*0 

250 

260 

35  0 


NT 


JF^INITE^EG.  1)  GO  TC  90C 
IFTIGAS. E0.2)  ITMAX  =1C 

fia*.  j 

start  of  balance  iccf 

ITER  =  ITER  4  1 
NHIN  =  FAIT  -  1 

CC  19 G  HE  =  1.  N0ANC 
LCCK  CUT  ENCTAU  AKC  TTF  ARE  ECUIVALE 
TTHL=SIHF<ENCTAUC 1,NB) ,1.CE45 ,1) 

CC  182  N  2  =  2  ,NWI N 
N 2P=N2-1 

TTHN2=SI NF (PNCT AU  (NZK,NE),1.0E45»2) 
TTHCNZM, NB) =TTHL 
TTHL=  TTH  N  2 

WhTJMIE  ,NE)=TTHL 

SUN=C.O 

IF  (NLTE. LE .3 1  GO  TC  186 
CC  18*  N2M.NMIN 
SLF=SLM4-TTP*NZ,NP) 

Ttdohmnei  =SUK 

TTUP(NB)  =  G  •  C 
CCWN C  NB) sSLF 
UFSTA  (NB)  =  C.O 
CCNTINU  £ 

CONTINUE 

NIAL  SHOILC  BE  A  ECU  EE  CF  TWC 
NT AL=E 

WRITE  (3,  1ICC) 

FCBKATT*  VIEFOP  *1 

<3,  6? "Hf  PPOF  <NZ,NL)  ,NL= 1 ,NLEVEL> 
WRITE <6,81  (VI8F0F  <12,NU  ,Nt=i,NLEVEl) 
WRITE  16.  8)  (VIBP0F(32,Nl)  ,NL=1  ,NtEVEU 
WRITE  <3,  1CG1) 


< K I T  E  <  3 ,  *  I  C  KP  l 
iRITE  (3.  1008) 
:CRM  A  T ( 1 F  6 1 12  < 
GO  TC  C  21f  ,21 


*) 


WKITET3, *) CRP(NP)  ,NE  =  1,WEANC1 

#  'IC«T'  •“t-’M.lM 

j 2 20 , 2 2 C  ,  2 A C  ,2 5 0 , 2 ( 0 , 260 )  ,  IGPS 

CONTINUE 
CO  TC  350 
CONTINUE 

CALL  CAL  Z  (fC2,NIAl, AC, \  IP) 

GO  TC  3  5  C 
CONTINUE 

CALL  CALZ (F2C,NI AL, ACtVIP) 

GC  TC  3  5  C 
CCNTINU  E 

CALL  CAL 2CFNO,NIAl,AC,VIE) 

GO  TC  35C 
CONTINUE 
GC  TC  3  5  C 
CONTINUE 

CA^Tpy  CCT.MAL.AC.VU) 

CONTINUE 

°W6?  iJft5.ll  N8,KC 

W  F  =  ‘IOC  (NEC,  100  ) 

KLN  =  KLCNEI 
KFN  =  KETNE) 

E\.N=BV  (1  ) 

bM  =  iJ^^^ir<l-^^^VE(NF,/29E.O) 

CC  360  N  2  =  1,  NALT 

ENOTAUf NZ,NE)  =  T A L F A X C KLN  , KFN , RP, P VN , TEMP ( N Z ) , BNCL  ,SNB « 

i6o *  c^ffssr^” 

'  ‘  CONTINUE 

IF  (ITFR  .IT.  ITFAX)  GC  TC  1P1 
CONTINUE 

IF  <NLT£  .LE.  2)  GC  TC  92C 
CC  91 C  N i  =  1,  NALT 

WRITE  (A, 8)  <VTBFCPtNZ,NL>  ,  NL  =  1,  NLEVEL) 

910  CCNTINUE 

CC  915  N  2  =  1  .NALT 

WRITE  (**,?)  lENTTAL  <N2,N'E)  ,NB=1  ,NPANO  ) 

CCNTINUE 


<0  0  0 
GOO 


915 

520 


VA. 


1C 


10 


30 


40 


C 

50 


C 

60 


ICC 


1 5  C 


STCP 

8  FCPMAT  (1F8E12.4) 
END 


SL0FOITI  fE  C8L7C GfSjKie^AC^VIE) 


TUS  SU8RCL  TINE  CALI 
CIN^NSICN  V  IE (4  0 , 2  ) 


FCP  F CP  A  GIVEN  GAS  AT  ALL  ALTITUDES 


/;^gjN.as5;is^sn;fU?^un!tE§?E?S8f^ssuu:!;ftf^fc!‘S5. 

3  iCCNC  (2  01)  ,CH£ML  (201) 


EXTERNAL  GAS 
WRITE  (6,  *)  NZ,NIAL,AC 

Hit  CALVIE (GAS,N2,NTAl,NEXT, Ar,VI«> 

NZ  =  2*MAL4N2 
N I AL  =  N  EX T 

IF((NZ42rNEXT)»GT  .NALT)  NIAL= (NALT-NZ) /2 
IF (NIAL. LE.  C)  PETLFN 
GC  TO  10 
END 


TH?SC^uSF^L^lkET?lL^LtA?is  VlEL?CPX?^R  fHR^  ALTICUDES  AND  INTERPCLA 
BETWEEN  TNE  END  POINTS  TC  CHECK  ACCURACY  AND  CHANGES  THE  AlTIUCE 
INTERVAL  IE  NFCESARRY 

NOTE*  yPSTfl.OCWN,  TTUF.7TCOWN  NU'T  BE  SET  BEFCRE  CALLING  THIS  *>CLTINE 
CIHENSION  VlBf40, 2), ATNL(IC) 

CCMMCN  VIEFCF( 20 1,25) 

CCNKCN  /C  EC  AT  A /  N ALT • N  LEVEL  ,NEANCtZFAO,KB(40)  «E  A (40) ,GUF(40) 

2  *C-DONN(  1 0 )  ,TTUP(A0*  ,T  ICC NN  (4  0  ) ,  UPST  A  (4 0  )  ,  COWN  (4 0  )  ,T  TK  2  Cl ,  40) 

3  »rONC(2Cl),CHEML(?Cl) 

EXTERNAL  GAS 

WRITE  (3,  ICC)  NZST jNIAL  .NEXT  .AC 

PCRNAK*  CALVIB  HERE  WITH  N2ST=*,I5,*  MAL=*,IE,*  NEXT=*,I5 
2  ,*  AC=*,F5.2) 

IF  (NZST. EC.l)  CALL  G AS  ( V ie ( 1 , 2 ) , N ZST  ) 

CC  10  NL  =1 ,NLEVEL 
VlB(NLti)=VIE(NL,  2) 

CALL  TTSET(NZST,N1AL,1) 

N  Z=NZST»  NI AL 

CALL  GAS  <VIE(1,2)  ,NZ) 

CALL  TTSET(N7,NIAl,2) 

NZ=NZ4NI AL 

CALL  GAS  (VIE (It  3) tNZ) 

CNFCK  INTERF  VALUES 
S  LN=  0  < 

CC  30  NL=ltNLEVEL 

ATNL(NL)=APS(ABS(VIE(NLtl)4VIE(NLt3)),.5-A6S(VIB(NLt2))) 

SLM  =  A TNL  (NL  >  +  SUM 
ATEST=SUN/(FLCAT  (NLEVED) 

WRITE  <3t  ISO)  ATEST 
FORMAT (*  ATEST=*,  F10 .3 ) 


If  pcssiele 


RT  =  FLCAT  (M  AL  )  / 1  • 

^^LC^KNT.tNE.RT)  GC  TC  5C 

CC  40  NL=ltNLEVEL 
VIE(Nt,3)=VIE(NIL, 2) 

CALL  TTSET(NZSTtNIALtl) 

NZ=NZST4NIAL 

CALL  GAS  (VIB(1,  2)  ,N'Z) 

GC  TC  £0 

CC  QUADRATIC  FIT 

STEP=NIAL*2 

STFP2=STEP*STEPP2. 

CC  60  NL  =  1,NLEVEL 

AFIT=  IVlE(NL,l)»VIE(NL,3)-2.*\IEtNLt2) ) /STEP  2 
BFIT=(VIE(NL13) -V  1 1 ( NL , 1 ) ) /STEC 
CFIT=VIE (NL ,2 ) 

XFI0=NZST4NIAL 

wtrmmhL 

X=NZ-XMI C 

THE  VIB  FCF  IN  VIE<I,J>  IS  LCG  CF  OENSITY 
VI6P0F (N Z , NL) *EXP  ( ( AFIT*X4EF I T ) * X *CF IT) 


5?SET_yPEJA,^OWNi 


CALL  TTSFT»N7ST.?*NTAl,11 


39 


CC  7 0  N8  =  1 » N  E  AN  C 
UFSTA <NB  )  =  TTL'F(NB) 

CCHN  (t\P)  =T7tCWN (  NE> 

msTsaftt  factor  /gains  cscilaticn  CP  2. 

IF  «  (A  TEST  *4.)  .IT./H  NEXT=NEXT*2 

RETURN 

ENT 

SLPPCITINE  TTSET  (NZST.NIAL.NCC) 

MCC  =  1  IS  TAl  FROM  L'FSTA  AND  CCWN 
MCr=2  TAL  FRCM  T  TLF  ANC  T  TTCH  N 

CCMMCN  /CFCATA/  N  AL  T.M  E  VEL  .N  EANC .  Z  F  AC  .KB  (4  0 )  ,E  A  (4  0  )  ,  CUF(  4  0  > 

2  ,r-OCNN(  4  C  )  »  TTUP  U  Cl  »T  TCCKN  <4  C  ) ,  LPST  A(4  01  » COWN  (401  «TTN  (201.40) 

3  ,CCNC (201)  ,CHEML  (201) 

IF  (MCC.NE .1)  GO  TC  20 
CC  10  NB  =  1  jKP ANC 
TTUP(NB) =LFSTA(NB) 

TTOOMN(NE»=CCNN(NE) 

N  2F=N  ZST 4NIAL-1 
CC  30  N?=N?ST.NZF 
CC  30  NB=l,NBANO 
TA=TTN (NZ,NE) 

TTUP  (NB)  =  TTLP (NB)  4 T A 
TTDOWMNf  )=TTCOWN  (NE)-IA 
RETURN 
END 


NLSCI'CTNETNUMBER’cFSfHESLEVEL,  NPS  IS  THE  NUMFER  CF  THE  BAND 
CINFNSICN  GNCSTE ( 2  C 1 > 

?CNMCN  /cFCA:A/0N^LflNLEVFL,NEANC,7FAC.<B(40)  ,E A ( 4 0 ) , GU F ( 4 0 ) 

2  ,  GDC  KNUD  )  t  TTUP  (4C)  .  IT  CCWN  <4  C)  tUFSTA  (4  0)  *  COWN  (4  0)  »TTW  (2  Cl  .4  0) 

3  ,CCNC (2 01) .CHEML  (201) 

EGUIVAL  ENCF  <G N C S TE ( 1 >  » V I6FCF f 1 , 1 ) > 

1PT(XFER( VIBFOP( 1, NLS) , TTHfl.NFSI  ,N7,  1, I, KB  (NFS  )> 

2  4  XFER  (V^EFCP(NZ,NLS  )  ,TTW(N7«NES>  .  NALT.NZ,  2, KB  (NES  )J  )*EA  (NBS 

3  4GDCWN  (NFS)  *S«1  (KR  (NPS)  .7FAC*TTC5CHNf  KRS»  ) 

4  ♦GUF(NES),’SN12(KE(NES)  .  ftUF  (NPS>> 


1PT(XFER( VIBFOP( 1, NLS) , TTHfl.NFSI  ,N7,  1, I, KB  (NES  )> 

2  +  XF£R  (V*EFCP<  NZ, NLS  >  iTTW(N7,NES)  , NALT ,N2, 2,KB (NES IJ  >*EA (NBS ) 

3  4GDCWN  (NFS)  *SM1  (KR  (NPS)  .7FAC*TTC5CHNf  KRS»  ) 

4  ♦GUF(NES),’SN12(KE(NES)  ,TTUF  (NPSM 
RETURN 

ENC 

SLFPCLTI NE  CC2 ( V I F  I  •  N2 ) 

§fP^I5?oH5E(«!JmSnUm!ao.).yvl;siG(?o1.,vVosic.?oi. 

2  ,VVNSIG  (2C1) .VVRSIG (2C1) 

PrKMfK  l  T  PFCP  (  ?  (]  f  o  C) 

3  ,CONC(2Cll ,CHEML (201) 


TVS  =  TVSIG(NZ) 

VTS  =  VTSIC-  ( N7  ) 

m  l 

FCP7  x  VIEF0P(N7,7) 

FCP2  =  ((TVS  ♦  RTFUN(N7. 2,  1) )*GNCSTE(NZ) 


)  4  RTFUN  (NZ  *  F»  ?)  I 


rvjr  _  _ . ,  .i  3i  i  .  .  .. 

!  :  ^:ifmi :  :  (in 

VIRI  (2)  x ALOG (POF2 ) 

FCP«  =  VIPFCF(NZ.e) 

FCP7  =  (UVS*0.2S  4  GLF  (  7)  *SN12  (KB(  7)  »TTUP(7)  ) 

1  4  GCOHN  !  7 ) *SM 1 ( KE (7) « TTCOKN (7) ) )  *  F  CP  2  . 

2  4  EM8  )  * FCP8  ♦  EA(15)*VI8FCF(N7*B)  4  E  A  ( 16  )  *  VIBFCP  (N7»  10)  ) 

FCP5  x  VIEFCP(NZ,51 
FCP6  =  V  IBFCP( N2  »€ ) 

FOP3  =  (FCF2M0. 2SMVS  4  GUP  (2) *SN1 2 (Ke (2 > ,TTLP ( 2 ) ) 

1  4  GCOHN  (2)*SH1!KEI2)  .TTCCHN  (2)))  ♦  PCP5MEM4)  ♦  V IS ) 

2  4  EA(10)*PCP»  4  EA  (13)*VIEFCF<N7.B  )  4EAC  14 ) *VI BFCf* <NZf  10) ) 

3  / ( E  A  ( 2 )  4  VTS  4  TVS) 

2  4  FCP6MEA(6)  ♦  VTS))/(EA(3)  4  VTS  4  TVS) 


:  (III 


♦  GLP<4)  *SP12<Ke<4)  ,TTUF<4)) 

1  ♦  GC0WM4)*SMl<Ke<4)  ,TTC0KM4))) 

l  »• 

VIEI<5>= ALCGCP0P5)  _ 

f  CP6  =  <FCF4M0.5*TVS  «  £UP ^ >?SN1 2 <KE <§) ,TTLF<6 ) ) 

VIBI (9) =  ALCC<<VIEFGF<Pi2,3)MVVNSIG<KZ)  4  *TFl M N Z,9  ,1 3) > 

1  ♦  GP0STE<NZ>*RTFL'MP2,9,11>  ♦  VIBPCP<NZ.7)*RTFUP<NZ,«,  15)) 

2  / (MRS I G  <PZ )  ♦  €  A ( 11 )  ♦  E  ?  C 13  I  ♦  EA<15»II  ... 

VIBI  (10  »=ALCG<  (VieFCP  <KZ,7)*  (VVNSIG  (N’7*  ♦  RTFIPI  ( N7,  10*16)  ) 

1  ♦  VIBFCF(KZ,3)*FTFUN<PZ,  10*14)  ♦  G  PCSTF  <  NZ ) *RTFUK  <  PZ,1 C» 12) ) 

2  / <  V VRS  IG < PZ )  4  €A  < 12  )  ♦  EA<14)  ♦  EM  16)  >  ) 

VIBI  (II  =  CCNC ( N  2 )  -  (FOPC  4F0P6  4  POP5  4  POP7  4  P0P4 

1  4  PCP3  4  FCF2) 

IF  ( VI B  I  <  1 )  ,LE.  0.0  VlPIfl)  =  0. 05*CCNC<NZ) 

VIPI1  l)  =  flLCGIVIBICl) > 

RETURP 

EPC 

SLERCIJTIPE  H2C<VTCI.PZ> 

CIFENSIOP  VieT<40l,€NCSTE  (2C1) 

CCPMCP  VIEPCP< 20  1,2!) 

CCPKCP  /  CCL  /  TVSIC  I2C1) ,VTSIC <201) ,VVUSIG<201) ,VVOSIG (201) 

2  ,VVNSIG <2C1),VVRSIG (2  (1) 

CCPMCP  /CPCATA/  P ALT , P IE VFL  , P P AND , ZF A C ,K8 (40 ) ,E A <40 > ,GL F < 4 0» 

2  ,GCCPN{ 40) ,TTUP<4C)  ,T  T  CCKN  <4  C )  ,  L'PST  A  <  4  0 )  ,  COM  <4  0  )  ,T  TP  <  2  Cl ,  40) 

3  ,CONC  (2Cl),r HE PL  (201) 

ECUI VAL  EPCE  (GNCST€(l)«VIBFCFfi,l) ) 

VIBI <2) =  ALCG  <  <<1VSIG<PZI  4VVKSIG(NZ) 

1  4  RTFUMP7,2,  1)  )*GNCSTE<NZ)  4VIEFC M KZ,3)*E A <5 > ) 

1  4  RTFUP<PZ,3,5)  )  *  V  IEFCF  (Pi Z  ,  2  )  4 

2  FTFLN< PZ,3,2) *GPCSTE <P Z)  ♦  VIPFCP < PZ,6) *EA <1 1) ) 

3  /  ( E  A  ( 2  )  4  EA<  5)  4  VTSIG(PiZ)  4  VVNSIG(NZ))) 

V IRI (4) =  ALCC<<RTFL'N<PZ.4,3)*CPCSTF<  PZ> 

1  4  RTFUP(PZ,4,6)*GNCSTE<P7))/<EA<3)  4  EA(6))) 

VIBI  (5) = ALCG<<RTFLP<P 2, 5,4) *  CPC STE<PZ) 

1  4  CTFLP  <P7,5,7)*VIBFCP<NZ,2))/<EA  (4)  4  EA<7)>) 

2  4  Ef(i;ni 


EA  <11  ) 


VIBI <7) =ALCG<(RTFLP(PZ,7,c)*VTPFOp(PZ,2) 

1  *  RTFUN  <P  2,7,13  MGNCSTE  <P7> )/ <EA<9)  4  EA<131)) 

VIST  <«)  sALCGMRTFL'N  <Ni,fi,10)  *VIEP0P  <NZ,2) 

1  4  RTFUP <N2,fl, 14) *GNCSTE <N7)  )/<EA< 1 C)  4  EA<14))) 

VIpI(l)=ALCG(CONC(PZ)-  <  VIEFC  MNZ  ,?)  4  VIBPCP  <NZ,7)  4  VIB  FC  F  <KZ,  6) 
Z^^IBFOF  <PZ,5)4VlEFCF<PZ,4)4VieFCF<N2,3)  4VIBFC  F (N  Z ,2 )  )  ) 

EPC 

SLBRCLTI PE  RNC<VIEI,PZ) 

C1PFPSICP  VIPI<40),GPC!TE <2011 
CCPMCK  VIEFCF<  20  1,25) 

3c555sii^sis(.:yvsiii?!H!y,?I6,2"i,-vvusi';,20,’-''''DSic'2'1’ 

CCPKCP  /CFT  A  TA /  PALT,K  LEVEL, PEA NC,7FAC,KB<40),EA(40),GUF(40) 

2  ,GDCfcN<4C)  ,TTUP<4C)  ,T T  CCWM4  C 1 ,  UPST  A <  4 0  >  , COKPI  <40  )  ,TTH  (2  Cl,  40) 

3  ,COPC<201) , CHE ML  <2C1) 

IfSiTSVfSEEc 

1  4  GLP<  1)*SN12 <KE<1) ,TTLF<1)  ) 

6  4  <  >FEF<VIPPOP  <1,2),TTM1,1  ),K2,1,  1,KE<1)  ) 

5  4  (EA<  3)  4  VTSIC(PZ) ) *  V IPFCF <N Z,3 )  4  2 . 0*CH EPL < NZ) ) 

€  /  <  E  A  <  1  )  4  VTSIG(PZ)  4  V  VFSIC  <N  2)  4  TVSIG(NZ))) 

VIBI <3)  =  ALCG <<  <TVSIG(PZ) 

1  4  GIP< 3)*SN12 CKE  (3) ,TTLF<3) ) 

2  4  GLOf 2)4SN1P<KE<2) ,TTLF<2)  ) 

2  :  SE3Htl<i*I8iIKi)1):nE8tin!U  )  *  VIEFCP  <NZ ,  2) 

5  4  2.04CPEPL<NZ)  )/<EA  (2)  4  EA<3)  4  TVSIG(PZ))) 

V IBI  (1 )  =  ALCG(CONC(P'Z)  -  V ieFCFCNZ,?)  -  VI BPOP  <N7 , 3)  ) 


oon  oo 


u 


RETURN 

£►0 


C 


n 


2 

2 

3 


1 

8 

3 

4 


SUBROUTINE 


_ C3(VI0I, NZ) 

C 1HENSI0 N  VIBK40  I  ,G 


CCHMCN 


GNCSTE1201) 

VIEFCP(201,2*) 

- 1201) fVTSIG (2C1) ,VVUSIG(201),VV0SIG(201) 

EVEL  ,NEANC,7FAC,Ke  (4  0  )  ,EA  (40  )  ,GUF(40) 
,GD CKN ( 4  C )  ,TTUP140)  ,7TCCWN  14 C)  ,UPST Af 4 C ) , CO WN (40) ,TTN 1 2 Cl, 40) 
CGNC (2  01 ) , CHE ML (201) 


CCMMCN  /CCL/  TVSIG  (201 )  »\ 


jiUNiUHiiintrL  *  £  u  I  f 


♦  GLP(  n*SN12(K€(l),T7UFin  ) 

♦  IXFEF1VIBPOP  (1, 2),  T  TP  11,1),  NZ.  1,1,  KB  (1)1 

4  >  FER  (VIEFCP(NZ,2),  1 7H (N? , 1 ) , NAL T  ,N 2,2, KB  1 1 )> ) *EA (1) 
4  GC0HM1)*SH1(  KC(  1)  ,  TT  CCKN  1 1 )  )  )* G  NCSTE  ( NZ) 
/(VTSIC(NZ)  ♦  E  A  ( 1 )  4  VVRSIC(KZ))) 


VIEI(  1)  =  ALCG(CONC (NZ) -VI EPCF (N 2,2) > 

RETURN 

ENC 


oo! 

€! 

I 

I 

€( 


FUNCTION  XFER  1 R A C , TTH ,N TOP ,N EC T ,K, KB) 
CCMMCN  VIBFCP1 2C1,2E> 

CIMENSICN  RAC(161)  ,TTK161) 

NINT  =  NTCF  -  NECT 

NFL  =  NINT  4  1 

IF  (NINT  .GT.  0)  GC  TC  100 

X  FER  =  C.O 

RETURN 

ICO  j;£NTINUjj 

TP  =  0 1  C 
EA3  =  0,C 

IF  (K  .GT.  1)  GC  TC  2  C  0 
NF  =  NTCF 

RFAC  =  VIEFCF(NFfll /RAC(NF) 

SUM  =  O.S'RA 
CC  1 2 C  NZ  =  1,  NINT 
NY  =  NINT  -  NZ  4  1 

TM  s  TTN  (NY) 

TP  =  TB  4  TM 
JE6l_  =  S  L  G3  (K E ,  T B  ) 

RE  =  RF AC*RAC(NY)/VIBFCF(NF,1> 

SUM  =  _SLM  •  (RB  -  RA)  *  (  E  B  3  -  EA3)/TM 


120 

ECO 


RA  =  PE 
EA3  =  E 


E3 


RE*SM2fKE,TE))/’RFAC 


CCNTINUE 
X FER  =  AES  (SUM 
RETURN 
CCNTINUE 
NF  =  NBCT 

RFAC  =  VIBFCP1NF  ,1) /R  AC(1) 

SUM  =  0  t$*RA 
Cf  22  C  N  *  =  2,  N'PL 
TM  =  TTMNZ  -  1) 

TB  =  TB  4  TM 

ttw  ^$3r’Te) 

RP  =  RF AC*RAC(N7) /VIEFCF  (NF,  1) 


SUM  =  . 
CA  =  Rl 
E  A  3  =  I 


iUM  4  (t»B  -  RA)MEe?  -  E A3)  /T N 


:20 


:  E3 

CCNTINUE 

XFER  =  AFS  (SUM  -  RB*SM2(KR,TE))/RFAC 

RETURN 

END 


FUNCTION  SL0  3 ( K , TAU) 

COMPUTE  CNE-HALF  TME  VALUE  CF  THF  FUNCTION  L  SIB  NOUGHT  THPEE. 

X  .EC,  1  ,  PARALLEL  E  AN  C 
K  .EC.  2,  FERPENO ICULAR  EANC 

T  =  ABS (TAU) 

IF  (T  .EC.  C.CI  GC  TC  300  42 

TP  nr  -ff.  nr  ir  ?rn 


-  .  ^ _ _  .  >'  V  \ 


22  0 


■30 


a 10 


1 50 


300 


IF  n  .GT.  i.C)  GC 
X  =  (T  -  G.5)*2.  C 
SL03  =  C CC2.583E-4*X 


TC  110 


6.014SE-7)*X  -  0 

1  return2 '3e??1*flLCG  m*T  4  (*<,C2S5C'?,*T 

110  CONTINUE 

IF  (T  . GT.  4.01  GC  TC  120 
X  =  IT  -  "  .5  )  /  1  •  5 

CL 03  =  -  3.51E£-4*X)*X 


,1027233)  *X 


-  1.  C464E-3)*X 


♦  0  .0  2195  C2  >  *  X 


1  ♦  2. 76  18 c- 3 ) *  X  -  7.66t29E-3)*X 

2  -  C.0626550)*X  4  0  .1 6  7C  7  0  2  )  *  T 
RETURN 

120  CONTINUE 

IF  (T  .GT.  10.0)  GC  TC  13C 
X  =  <T  -  7 . 3  )  /  3  ,  C 

SL03  =  C <  C  ((  <1.  7lccE_t  «x  -  4  .99S1E-4) *X  4  0  .  C  Cl  1  377  2)  *X 

1  -  C.00 334638I *X  4  0.C1GC98£1)*X  -  C. 03 090774 )*X 

2  ♦  C. 10 *86302) *T 
RETURN 

130  CONTINUE 

Y  =  ftccm 

IF  (T  .GT.  100.0)  GC  TC  140 
X  =  S .  2  1C34C371976/Y  -  ?.  C 
SL03  =  <<(<(?. ETE-5  -  S.78E-5*X)*X 

1  -  0  .  CO  1 285  8) *X  4  C.0340267) *X  4  0, 

RETURN 

140  CONTINUE 

IF  CT  .GT.  ICOO.C)  GC  TC  15C 
X  =  2?. 6310211 1 56  3 / Y  -  5.C 
SL03  =  ((  CC8.3E-6  -  2 . 01E-5*  X) *  X 
1  ♦  C.C1263C6)*X  4  C.2726344)*Y 

RETURN 

150  CONTINUE 

X  =  13. fl 5  51 055  766 /Y  -  l.C 

SL  0  3  =  CCC  t  <C5.36€E-4*X  -  C.CC15638)*X  ♦  O.OC16306)*X 

1  -  0.0gC1371)*X  -  C.Cf  2-*757)  *X  4  0  .  C3  01  34  1)  *  X 

2  ♦  C.2315?20)*Y 
RETURN 

2  C  0  CONTINUE 

IF  (T  . GT.1.0)  GC  TC  210 
X  =  (T  -  C.E)*2.C 

SL03  =  (If  1. €  07  E*4*X  -  i.0S6EF-3>*X  -  0.08806C5)*X 
1  ♦  C .09  74C7381*  ALCG  CT  )  *T  4  0 .4161293) *T 

RETURN 

210  CONTINUE 

IF  CT  .GT.  4.0  GC  TC  220 
X  =  CT  -  2.51/1.5 

6103  =  (<(<<(( 5 .6C4E-4  -  2 . 764E-4*X ) *  X  -  8.  066E-4>*X 


4  2.45  8E-4)*X 
3210431)  *Y 


-  1. P43E-4I *X 


5C4E-3) «X 
!  3  C7 133 ) *  T 


4  0. C197235) *X 


GC  TC  230 


4  2.22*4E-3)*X  -  6.1 
f.-TCj06  3  1E€2)»X  4  O.i 

CONTINUE 
IF  (T  .CT.  10.0) 

X  =  (T  -  7.01/3. C 
SL  03  =  (( (  ((  Cl  .476E-4*X  -  4.261E-4MX  ♦  1.0472E-3)*X 

-  3.31 C£E-3|*x  ♦  0.01C7584) *  X  -  0.t3c6259>*X 
_ 4  C  .13 1 C340  ) *  T 
RETURN 
CONTINUE 
Y  =  /SLCCCT) 

IF  CT  .GT.  100.0)  GC  TC  240 
X  =  6.2  1 G  34  0  37 1 6  76/  Y  -  3.0 

EL  03  =  C(CCC(4.1fSE-4  -  C2.131E-4*X  4  3.79E-5  )*X) *X 
♦  7  .79  6-5 )  *  X  *  6.04E-5MX  -  O.C022C66)*X 
4  0  .03  30669) *X  4  0.4167794)*Y 
RETURN 
CONTINUE 

IF  CT  .GT.  1000. C)  GC  TC  25C 
X  =  27. 631C2111563/Y  -  5.0 

SL  03  =  C  C  ( 1.6KE  -5*X  -  3.642E-4MX  4  0.0134767)*X 
4  0 .37CCS86>*Y 
RETURN 


♦  9.4516*4) *X  - 

♦  C .32  54546  )*  Y 


4  12E-4* X ) * X  -  8.441E-4»*X 
3.47£6E-3)*X  4  0.0335987)*X 


c^nVinu  E 

SL 0  3  =  C.O 

RETURN 

END 


1*3 


A 


?USRH8MM2  <K’T<SU 

IF  (T  . EC.  C.C)  CC  TC  ICO 
IF  (K  .  (T.  1)  GC  TC  ZCC 
IF  <T  .  CT.  l.P)  GC  TC  lie 

SN12d=*  mm  wa  4E-1E*T  -  1 .  in7F-  1?)  *T  ♦  2.11117E-11)  *T 

-  €.43  16C9E-9) M  ♦  1 . 2 t 7£ ?E - e ) * T  -  7 .576  26E- 4) *  T 

-  0. 08  264246  4  C. C19P45€4*AICC (T>  )*T  ♦  1.0)*C.5 
RFTU'RN 

CONTINUE 

IF  (T  . CT .  4  . C)  GC  TC  12C 
X  =  <  T  -  !>.*) /i.5 

SKI 2  =  f U tt U 3. C0591F-3  -  1  .  41 4€6E - 3«X I  *  X  -  2 . 3  ft 08 5 F - 3  )* X 
4  9  .40  736E-3) **  -  C. C2’54C7)*X  4  0  .05230**) *  X 

-  0.1Gl5ce4)*X  4  0 . 1  ?  €2  743)  *  C  .  5 
FFTUFN 

ccnttnu  E 

IF  (T  . CT.  1C. 0)  GC  TC  13C 
X  =  (T  -  '».P> /7.  C 

ENl 2  =  t  ( ( t< <(5  .  C8470E-4  -  1  .90  1  €1 E-4*X) *X  -  E. 27^9  4F -4 ) *X 
4  2 .24  274E-3) *X  -  5 . 5  4  2  S  8  F  -  3  >  *  X  4  C . 0 1 29 Cft 0 4  )  *  X 

-  0.02979C75)*X  4  C.C£56?2F?)*G.F 
FFTUFN 

CCNTINU  E 

IF  (T.GT.  100.0)  GC  TC  14G 
Y  =  ELCC(T) 

X  =  9.21C34037197F/Y  -  3.C 

EN12  =  (<((((  7.?7€E-£.«X  -  1  .  El  81E- 2)  *X  4  9 . 0 12 8E-4 ) *  X 

-  6.49fE-4)*X  -  l.C€7ZE-3)*X  4  0.  0  14936M*X 

4  0.43  iP273)*G.5/T 

RETURN 

X  =  9. 2  1  G  3  4  C  371 9  7€ / Y  -  1  .  C 

EM2  =  <m5.C74€£-2  -  4  .  €  8 1  E-4*  X>  *  X  -  1.584PE-3)*X 

4  0.01  1752F)*X  «  0.4C79C74) *0.5/T 
RETURN 
CCNTINUE 

IF  (T  . GT .  l.C)  GC  TC  210 
3=2*. C729185*T 

SN1?  =  <<<  <<t(2.71E-l€*T  -  6  *  P3  72E-  14  )  *  T  ♦  1  .  5 8 3 1  4E - 1 1  >  *T 

-  3.3qe73E-9)*T  i  7 . C  1  € 1 8 E - f ) * T  -  1 . f 1 59 4E - 4 ) * T 

-  0.06c£453  4  0.014942f4*eLCG(T) »*T  4  1.0M0.5 

RETURN 

CCNTINUE 

IF  <T  .C-T.  4.0)  GC  TC  220 
X  =  <T  -  2.FJ/1.E 

SN12  =  ( ( (  t  (  ((?.  22E-3  -  1,  12  C3E-3*X  )*X  -  3.04P9E-7MX 

4  7.93  45E-7)*x  -  b.Q213281)’X  4  0.C521C86)*X 

-  G.1213EEC)*X  4  C.2ECP902) *0.F 
RETURN 


CCNTTNU  E 

IF  (t  . CT,  1C.0)  GC  TC  23  C 
X  =  (T  -  7.01/3.  C 


SN12  =  <  ( <  ( < <4  ,727E-4*X  -  1.2455E-3MX  4  2.  52E3E- 3)  *X 

-  6.73?3E-3)*X  4  C.C174769)*X  -  C.C421  447MX 

4  0.C95E415)*0.5 
RETURN 
CCNTINU  E 

IF  (T  .GT,  1C0.0)  GC  TC  24C 
Y  =  «LOC(T) 

X  =  9.2  1C34G371  97F/Y  -  3.0 

EN1?  =  i(((9.1F-c*x  -  1.909€E-3)*X  -  1.8059F-7)*X 
♦  0.12F?129)4X  4  C.£35?594)*C.5/T 
RFTURN 
CCNTINUE 

y  =  fiocm 

X  =  9.2 1C?407797(/y  -  1.0 

E N 1 2  -  ( ( <(€.444E-4*X  4  2.1  53F-4)*X  -  5.099E-4)*X 
4  O.021€1CO>*X  4  C.55105?I*C.5/T 

FFTURN 
CCNTINUE 
E  N 1  2  =  t  .5 
RETURN 
ENC 


FUNCTION  1 AUMAX  tKL»KE,RN,BV»T»FLAM,S»DEN> 

0*SCFT  jRN'/T)*FLAM*DEN*S*SJMAX<KL,  *e,SIG) 

RETURN 

END 


FUNCTION  SJMAX  ( KL.KE ,SIG) 


ini  n»MfKJf&p«G,H 0F  7HE  ;IR0MGEJ,  llKE 

GC  TC  (  110,120, 130  ,  KE 
110  CCNTINUE 

IJ2  5Jr°/JIC)  *  1-2 3 * 5 * * * * *’ 

GO  TC  2CC 
120  CCNTINUE 

7JM  =  0  .5MSCRT  t2.0/STG)  -  1.0 
A  =  I JM  «  l.C 
GO  TC  2  0  C 
130  CCNTINUE 

SJMAX  =  C.C 
RETURN 

200  CCNTINUE 

S  JM  =  A  *EXF(-SIG*TJM*  ITJN  ■»  1.0)) 

?JMAX  =  SJF/CJPAFTIKL.SIG) 

RETURN 
END 

FCjplR?fCOPUT€5  CTfiP^CTATICNAL  PARTITION  FUNCTION  FCR  A 
RIGID  PCTATCP,  SIG  .IE.  0.2. 


K  =  1,  ALL  LEVELS  FCFUL A  TEC 
K  =  2,  EVEN  LEVELS  FCPL’LA  TEC 
K  =  3,  CCC  LEVELS  FCFLLA1EC 


IF  (SIG  .r-r.  0.2)  STCF 
X  =  SIG 

IF  <K  -  2)  110.12C.1TC 

COEFFICIENTS  ARE  FRCN  RUN  CALTE-EV 
110  CCNTINUE 

CJPAET  =  U((O.0C28€ECO2?2»57*SIG  4  0 .0 12 S7l 8 E35 471 > *  X 

1  ♦  0.06  €€  756  1€4520|*X  ♦  C.  233332133  270  *X  4-  1 .0  0  COO  CO  CO  17  ) /SI  G 
RETURN 

120  CCNTINUE 

2  ♦  f.02S737Ciq36Ae)*SIC  ♦  0 . C 32 8350 609506 ) * STG 

2  4  0.16  €670220  770) *SIC  4  0 . 4 S E« 5999 94 56 > /SI G 

F FT URN 

130  CCNTINUE 

CjPAFT  =  < (t << M 29.57C211707E  -  62 . S 8 7459E5 02*SI G ) * SIG 

1  -  5.495763736651*SIG  4  C. 505 75C942 701) *SIG 

2  -  C.C1 7C2S7E51£7ei *SIG  ♦  C»C32C315C8C069)*SIG 

3  j.*T0j,16€6E311238  7)*SIG  4  .  5  C  C  C  CC  00  0  f  4  4)  /S IG 

ENC 


FUNCTION  SIMP  <A,F,K) 

LlMENSICN  A  ( 3) 

IF  (K  -  2)  IOC,  110, 12C 
CCNTINUE 

SIMP  =  18 . 0* A ( 2  >  45.0'A<1>  -  At3))*H/12.0 

RETURN 

CCNTINUE 

SIMF  =  (E.0*A(7)  ♦  5. C* A  (3 )  -  A(1))»H/12.0 

RETURN 

CCNTINUE 

SIMP  =  t A  •  0* A ( 2 )  4  A ( 1 )  ♦  At  2) ) *H/3 ,0 

RETURN 

END 


>  v  \ 


110 


120 


130 


IF  <7  •  EC  •  C.O  GC  TC  360 
IF  (K  .C-E.  2)  GC  TC  200 

iw^.sin.i00  ic  110 

SMI  =  C  M(12.354?CE-C3*X  -  8  .50535E-03)*X  *  2 
1  -  7 • 07 { 77£- 02)  *  X  •  l.E3735E“Cl)*X  -  3.01749E 

RETURN 

5FnIInVIt.  10..,  CO  ,0  1?C 
X  *  IT  -  7.0/3. C 

SMI  =  M(M5.q22«7E-C4«X  -  1 .79043E-03)  *X  *  4 
1  FETU*S<*n°E"C2*  **  4  2.*91CSE-t2)*X  -  E.0170CE 

CONTINUE 

IF  IT  .GT.  100.01  GC  TC  130 
Y  =  ALO  G  ( T ) 

X  =  5.2 1034037197E /Y  -  3.0 
SMI  =  CC  <«<<-€. 24516E-C4*X  -  f  .75378E-  (4)  *X  ♦ 

1  -1.2815CE-03)  *X  -  E.31784E-C3)*X  ♦  3.78904E- 

2  FM*r*,E*,1,/] 

CGNTINUE 

X  =  l^9^?]io37ie7€/v  -  1.0 

SMI  =  1  Ml  1(1.  l2C38E-04*X  -  2.9  3474E-  (4)  *X  - 
1  4  1.30C37E-C2)  *X  -  3  .  f  0  08  8  E  -  C4 >  *X  ♦  3.234U1E 

♦  8.174G7E-C1) /7 


.56223 E-  C2)*X 
•Oll’X  *  0.426845 


.051  09E-  C3  )*X 
-02 l*X  4  8.13 


134123 


:oo 


:ig 


220 


•30 


300 


c^nVinue 

IF  (T  .GT.  4.0)  GC  TC  210 

X  =  (T  -  2 . 0  )*  0  •  f 

SMI  =  f  <(<<1.396?  IE-03  *X  -  5  .  2590 4 E- C7 > *  X  ♦ 

-  4. 740  35E-C2)  *X  ♦  1. 2 f 1 7rE- Cl ) *X  -  2.86333E- 
RETURN 

Wn\]W.h.  10.0)  GC  TC  22  C 
X  =  <T  -  7.0)/3.C 

SMI  =  <(<( f3.7642CE-04*X  -  1 . 253P4E -03) *X  4  3 
-  1.05360E-02) *X  4  3 . C 7476E - C2 ) *X  -  P.14705E 
RETURN 
CONTINUE 

IF  IT  .GT.  100.0)  GC  TC  230 

Y  =  ALOC-(T) 

X  =  5.2 1G340371976/Y  -  3.0 
SMI  =  <  H-7.53415E-C3*X  -  8  .  <»7l2*«E-  C3>  *X 
4  1.30t24F)/T 
RETURN 
CONTINUE 

Y  =  B  LOGO) 

X  =  5.2  1 C340  371 5  7F / Y  -  1.0 

'S12E-04*X)*X 


2.  98590E-0  3)  *X 
C2)  *X 


4.36163E-04)  *X 
-0  2)  *X 


1.  5  8245E-02)  *X 
Cl)  *X  4  0.517573 


.411E8E-  (3  )*  X 
-0?)*X  4  0.19413C8 


♦  6.  E2199E -02)  *X 


SMI  =  << <t(l.«9p67E-C3  -  *,8251° 
4  1  .21174F-(J’)*>  4  4.7EU5F-C2 
RETURN 


>*X  4  1.10445 


59665E-03) *X 


*)'?66 


CONTINUE 
SMI  =  l.C 


fEJURN 


1*6 


FROG  F  AM  E4IKT  (  CITFU7»TAPE7  ,  T  A  R£4 »  T  AFE6 ) 

f  issfir.fs  tt5T  F*tT  ir  ce5g,s  ecnc 

CIKFWSTCN  XIC  (4)  ,FAC  (26  >  .RLAM25) 

CIMFWSICW  l SCI? 5)  ,  ICC (26)  ,eC  M25  )  ,  E!V  (25 ) 

CIKEKSTC  W  2  AN G(71  ,  ST AU(7),'fiAC<8> 

Cl  KENTON  l  EC  (2  5)  ,L e C  (  2  E  >  ,?TR  (25) 

CCKMC*  / FCF CCM/  V  IEPCF  (  2 C  1 ,  1  C ) 

CCKMCN  T  ft'  (201,  25)  ,R  AC  <  2  CD  ,T  (201) 

CATA7  2ANG  IT) ,1=1,7) /0. C,20. 0,4 0.0,60  .0.70. C, 8 C.0, 90. C/ 

2 € 0  FCRMAT  ( 2 1 6 , 1 F5 E 1 2 . 4 ,2 1  2 ) 

26?  FCRMAT  (21€,1PE12.4,2I2,1X,4?6) 

R£AD(4,3  78>  7 

READ  (4,3621  TGAS.M7E 

RE  Arm, 3  1C)  ( XI C  t  K  1  ,  K  =  5  ,  4 )  ,NlfVEL,NE?ANC,AMASS 
3 1 C  FCRMAT(4Ae,ex,I2,2X,I2,3X,Fe.2> 

WRITE  (2,37?)  T 
37  r  FCRMAT  (1F6E12.4  ) 

WRITE  (2.36C1  IGAS 

WRITE  (  2,3621  NtEVEL»NEAND,ANASc,IGAS»NLTr,(XIO(K)»K=l,4) 

OC  381  Ij  =  1,  NIEVEL 

281  CCNTINUE 

282  CCNTINUE 

RE AO (4, 360)  IALT.-ALT 
WRITE  (2.260)  IALT,JAL7 
ft  FCRMAT(1F8E1C,4) 

118  FCRMAT  (1F8F12.4  ) 

DC  331  N2  =  1,IA1  T 

331  cPSIe’nI-i  m?PCMNZ,*U  tM  =  l,NL£VEL) 

336  WRITE  (7,6)  (V IOP C F (N 2 , N L 1  ,NL= 1  ,NL EVE L > 

CC  337  N2=1,IALT 

377  REAO(4,118)  (TAU(N7»NS) »N6=1»N6ANC) 

jS  FT  =  f> £A  N C 
CC  84 C  J=1,*SET 
NfC=LEC( J> 

Nl=MOC(NEC/lflO, IOC) 

CC  6  0  E  N2=1,IALT 

600  RAC(N21  =  FACCJ)*VIEFOF(f2,M  ) 

700  FCRMA7(1F1,4A6, 8HEANC  N C . , I  7 , 6X . 7PL A NC A  =  . 56 . 2  ,7H MTCRC  M,4X,4A6//1 
WEITE(4,7CC1(XIC(K),K=1j4),J,RLAM(j)  »(XIC»K)  ,  K  =1 , 4  ) 

IF  (NLTE  .CT.  11  CC  7C  707 

’04  FCpM  AT  (4  EH  ASSUMES  LOCAL  THR* KCC YNAM IC  ECU  IL  IERI UM/I 

705  WRITE (6, 7  C4) 

CC  TO  71  * 

706  FCRMAT  (28H  LACK  CF  THERMAL  ECUlLI°RItM  ,16/) 

707  WRITE  (6  ,706)  NLTE 

710  FCRMATd  1 8F  EXOATKCSFH  ERIC  VIEWING  S?  -  --  --  --  --  -  -ENC 

1CATMCSPHERIC  VIEWING---- - -  -  ) 

715  WRITE  (6,710) 

720  FCPM AT ( 1  1  6  W  TAN  HI  LIKE  BAND  RAC  ff  ALTC  ------  PARC  RADI 

1ANCF  AS  A  FUNCTION  CF  2ENITF  ANGLF  (WATTS/ CM 2- ETERI-  -----  ) 

WRITE  (6  ,720) 

720  FCRMAT(118W  (KM)  JKATTS/CM2-S7FR)  ft  (KM)  C  DEG  -  L’F  ?0  C EG 
1  40  CEG  6C  CE£  7C  CFG  80  OFC  60  CFG-HCRZ  F 

WRITE  (6,7’C) 

740  FCRMAT  (F5 , 0,1  PE  14 . 4 , 0 FF i ?. o , 1 X ,1 P7 El ? .4 ) 

KTAU  =  WCC(LGU( J) , 1C) 

KEET=IALT-1 
CC  830  K  =1, KSET ,5 
IF (K. CT. 61)  GC  TO  750 
H  T  =K  *59 
GC  TO  76 C 

750  HT=15C*( K-91) *5 
760  C1=HT«6371. 

C  T ► P , 6  ,IN  THE  CO  LCCR  MAKES  IT  CNLY  CALCULATE  VALUES  FCR  0  AND  «C  CE 

CC  81 C  L  =1,7,6 
M  1  =  K*  1 

STAU(ll=?.0 
ELIA  =  C.O 
$RAC(l)=C.C 

stmiih  * 


A  v  \ 


PIS 


PI  € 
PI  7 


CC  79C  K  =  K1  »HSFT 
IF  (N.CT. 51)  GC  TO  77C 

Ens'nr1 

770  C2=15C4( H-91)  *‘>♦€771 
780  AKG  =  ASIN(STMPHI) * C 1 /C  2  ) 

CS=SCm  C 1 *C 140 2* C2- 2.* Cl *C2* CCS (PHI -A PG)  ) 

C£LS=(0S-RAKG)*l.CE4t5 

RANG=CS 

CELTAU  =  CEL'M  0  .?*  (TAL(H, 

STAU(l)  =  s1SU(LJ  4  CELTAU 
SL1E  =  SL1(*TAU,S7AI(L)) 

RADA  =  SAC  (H-i)/TAU(H-l,J) 

FACE  =  FAC  (n/TAl  (P,J  ) 

'FAC  =  (RACE  -  RACAI/CELTAU 

Sc  AD  (LI  =  SR  AO  (  L  )  4  SFAC*(SL1A  -  SLIP) 

S  L 1  A  =  SL1E 
SAVTAl=STAl (L) 

790  CCNT IPUE 

SFAO(l)  =  SR AO ( L  )  ♦  RACE’SLC  (KTAU.SAVTAU) 
€10  CCNT  IMJE 
S  A VT  AL  =  0  .C 
STAU(7)=C.C 
SL 1  A  =  C.C 
VFANG  =  C.C 
SFAD( e>  =  c.o 
K'ET=TALT 
Cf  81?  M=K.KSET 
L  =  PS  E  T-M  4K 

IFIL.GT.ei>  GC  TO  Pll 

EPVS-S1I”1 

Pll  C2=1EC4( L-S1>*54* 371 
P 12  CS=SQFT( C2*C2-D1*C1) 

PSET)  GO  TC  813 

GC  TO  81? 

PU  CELS=(RAPG-C?)*1.  CE4  C5 

mu*  * C£LS 

-  =  S  T  A  U  (  7  I 

=  CELS*  f  0.?*  <TAlfL,jMTAt(L*l.J>>> 


SAV^fu 

CEL  T  AU  .  . . 

lHii,:,sCj««KiTai;«vilu 

FAOB  =  F AO  (L )/ TAL ( L ♦ J ) 

RADA  =  . . . 

SFAC  = 

'RAC  C8> 

?L1A  = 

CCNTIPUE 
P1=K41 
R  A  NG  =  G • 0 
KSET  =  IALT 
CC  8 1 P  HsKI.vset 
IF (M. GT »  S 1 )  GC  TC 
C2  =  M4  59* {371 
GC  TC  81  1 

C2=1?C4( P-91) *546371 
CS=SQF  T(  C2*C2-01*  Cl) 

C|LS=I0S-RANG>*1. FF4 C5 
VRANG  =  VRANG  4  CfcLS 
RANG=CS 

=  ST  AL’(  7  ) 

=  Cels*  < o  .?*  etauk,  j>  «taick-i»  j>  >  > 

=  ST  AIM  7  )  4  DELTAU 
SU  <KTAU,STAL<7I) 
RADJP-ll/TAKH-ltJ) 

FAP<P)/TAUf',j) 

(RACE  -RACA) /DFITAU 
=  SR  AD  ( 8  )  4  SFACMSL1A  -  SLID) 


RAC(l4l)/lAL‘  <L4l,j) 

(RACE  -  RACA)/rELTAU 
=  SR  AD ( P )  4  SFAC*(SLlf  -  SLIP) 
SL1E 


P 16 


SAVTAU 
CFLTAU 
S  TAU  (7) 
SL1B  = 
RADA  = 
RAOP  = 
SFAC  = 
SRAC  (8) 


P 1 8 


4U4 

(  30 
(40 


C?HlNUESL,E 

SFAC  (8)  =  SRA0(8)  4  R A CE*SL 0 (KT AU , S T AU ( 7) 
WRITE (3  »444A )  HT  ,SR AC  (  P ) , SR A C ( 1 ) .SR  AD (7) 

-  -  IP  3c 1 7.4 ) 

»3fAC  (8  )  «  HT  f (SRAC(M)  .  M»  1 


MW"'* 


CCNTINUE 
STOP  77 
END 


7) 


*  » \ 


FUNCTTCN  SLO  (  K » 1 A  L ) 

T  =  AES  (TAD 
IF  (K  .EG.  2)  GC  TC  2C0 
IF  (T  .GT.  4.0)  GC  TC  '110 
X  =  (T  -  2  •  0  )  *  0  •  5 

SLC  =  ((((<( C. OOt4110*X  -  C.Ctl  4  0.  0  0  683  7  E»  *X  ,,,  „  M, 

1  -  0.02<e0€C>*X  *  0.0752553)**  -  0.  D01  C?3)  *X  4  0.6S69377)*T 

RETURN 

110  CONTINUE 

IF  CT  .CT.  1C.0)  GC  TC  1?C 
X  =  (T  -  7.01/3.  C 

SLO  =  (  U(((0.00C151E*X  -  C.CC07254)*X  4  G.  0  C  230  7  6)  *X 
1  -  0  .Cfl8C248)*X  4  C.C274665DX  -  3 .  C9  1 2  08  2)  *  X  4  0.3*39255)*T 

FETURN 

120  CONTINUE 

Y  =  ALOC(T) 

IF  <T  . CT.  ICO. Cl  GC  TC  170 
X  =  S.21C34Q371976/Y  -  3.0 

SLO  =  (  ((< <<C. OOTCOCS'X  -  0.6C01530)*X  -  C. 0 0  Cl  4 9 C) *X 
1  4  G. 00  13543) *X  -  0.0CEIC*C1*X  4  (I.  C837C6?)*X  4  1.141*564) *Y 

RETURN 

130  CONTINUE 

X  =  5.2  103*0371976/3  -  1.0 

SLC  =  ((((((((0.(6714'  -  r.0fE222*X)  4  0.005211)*X 

1  -  0  .0  0  4  85  6 ) *X  -  C. 00*477) *X  4  C.0C7743)*X  -  C.01*9C0)*X 

2  4  0.122526) *X  4  C.54C111DY 
RETURN 

200  CONTINUE 

IF  (T  . GT.  4.0)  GC  TC  210 
x  =  (T  -  ?.0)*n.* 

SLO  =  (  ((( (<0.00C24  31  *X  -  C.C(1C5°1)*X  4  0 . C C 41  8  3 0)  *X 
1  -  0 .01  6076D*X  4  0.0*61765  )«X  -  0  .  20  1327CDX  4  0.716P944DT 

RETURN 

210  CONTINUE 

IF  (T  .CT.  1 C. 0 1  GC  7C  ?2C 
X  =  (T  -  7.0/3. 0 

SLO  =  ((((((  0.  0  0  0  11  (9* X  -  C.C00  4699  )*X  4  C.  00  16637)  «X 
1  -  0  .00  65  3  45  ) *X  4  0 . 0 2 *55* 2 > *X  -  O.10C3592)*X  4  0.42e2944)*T 

RETURN 

220  CCNTINUE 

lJ  .  o, 

1  -  0,05  68*5) *X  -  C  .  04  6  61  6>  *X  4  t.0  45072)*X  4 

2  -  O.GOC8*7)*X  -  C  .  G6  2C5C ) *X  4  C.163?1°>*X  4 

ENCU  N 


C6  2211 1  *  X 
C.  00C249)*X 
1.  3  997  65  )  *  Y 


C 

r 

c 

r 


no 


FUNCTION  SL1  (K.TAU) 

IF  ( K  .EC.  1)  CONFUTE  FUNCTICN  VALUE  FCF  PARALLEL  BANC. 

IF  (K  .EC.  2)  CONFUTE  FLNCTTCN  VALUE  *CR  FE R F EN 3 ICU L A F  EANC. 

Jf",«R?«!1$,  GC  TC  2  C  0 
IF  (T  .  C T  .1.0)  GC  TC  11C 

X  =  (T  -  0.5)*?.  C 

SL1  =  (J((1.1(’E-*’X  -  1.5«.*26-4>*X 
l  -  C.C3  **3725) *X  4  C. *  €  1 354 7  1  ) * T*T 
RETURN 
CCNTINUE 

I*  (T  .GT.  4.0)  GC  TC  120 
X  =  (T  -  2 .5  )  /  1 , 5 

SL1  =  (  ((((3.497E-H  -  5.4cE*5*X)»x 
C.0114G87)*X  -  0.06  1  3*64  )  *  X  4  1} , 


1  4  [  .0 

RFTURN 
120  CCNTINUE 


4  C.  00  2761  49)  *X 


-  7 , C3  28E- 3) *  X 
3531337) *T*T 


1 


IF  (T  .GT.  10.0)  GC  TC  DC 
X  =  (T  -  7,0/3. C 

SL1  =  (<(( (7.771 E-4  -  1.87EF-4*X)*X  -  2.9392E-3)*X 


4  0.0120C34)*X-  <$.0*Cil3*)*v  4 
RFTURN 

130  CCNTINUE 

IF  (T  .GT.  100.0)  GC  TC  140 
Y  =  AUOG(T) 

X  =  5.2  1  0*40  371  57E/Y  -  3.0 

SL1  =  (  (  (( (6 .67R-* *X  -  6.76FE-4)*Y  _ 

1  F4TCj,00  1  184  1)*X  •  C.CC56382)*X  4  C  .  6538  *6  5)  *  T*Y 

140  CCNTINUE 

rsiiift^cf),ooB,n  Gc  Tc  i5c 


0  .  235431  C)  * T *T 


4  0.0  0  21  934)*X 


49 


VS 


1 

150 


1 

2 

200 


1 


X  =  2?l  ES162111553/Y  - 
511  =  (  <<i.£?6E-4*X  - 
0 ,861€2<:e)*T*Y 


♦  0.8611 
RETURN 
CONTINUE 


5.  C 

1 .e6C4E-3>  *x 


1.6852E-3)  *X 


y  =  fitocm 

X  =  IX.  61551055  756 /Y  -  1.0 

511  =  (  l<<  U(2.56F-4*X  ♦  2.3C9E-3)*X  -  7.6  00£-3>*X 

♦  6  .5175*31 *X  ♦  2.37F-4)*X  -  l.?795E-2)*X 

♦  0  .030224) *X  4  C.P3P671) *7<Y 

RFTURN 

CCNT INU E 

IF  (T  .  CT.  1 .0)  CC  TC  210 
X  =  (T  -  C.5>*2.  C 
SL1  =  (  CU7.C9E-6 
-  0  .02  632252)  *X 
R PTUFN 
CONTINUE 

IF  (]  .  CT,  4, C)  CC  TO  220 
X  =  (T  -  *.c ) /I  .  5 

511  =  ((<<(’. 14E-4  -  3  .26F-5*)) *  X  -  1 .  3 29 CE - 3 ) *  X 

♦  «.3415E-3)*X  -  0.0530254) *>  4  0  .  78 01 694) * T *T 
RETURN 


*X  -  1,212!E-4)*X  ♦  1. 87580E-3) *X 
4  G.i£*£7333)*T*T 


220 


1 


240 


1 


250 


1 

2 


CONTINUE 

IF  CT  .CT,  10.0)  GO  TC  230 
X  =  (T  -  7.0/3.  C 

5L1  =  (  I C  C  (5 .23  0  £-4  -  1 . 1 c  2E  -  4  *  X  )  *  X  -  2 . 2 \">  0  E -3 )  *  X 

♦  0.01C2C27)*X  -  0  .0i542r2)*y  4  C . 271 8 1 0 1 ) * T * T 
EFTURN 

CCNT  INU £ 

IF  {  T  .  CT.  100  .  C  )  C-C  TC  240 
Y  =  £LOCm 

X  =  5.2  1C‘,403?197€/Y  -  3.0 

5L1  =  C  (  (  C (1 .06 P£-4*X  -  1  .  3224E-3)*X  4  7.P304£-J)*X 

♦  8.184Cf-7)*X  -  e.£ C 62GF-2  )  *X  4  1  .  C 57 7 8  1 )  * T * Y 


Gr  TC  25  C 
5.0 

i  2C1 2E-3) *y 
1715) *T*Y 


F  ETU  CN 
CCNTINU E 

IF  (T  . CT.  1  COO. C) 

Y  =  «LOC(T> 

X  =  27.  £310211  ltr53/Y 
511  =  C  ((2 ,€68E-4*X  -  c. , 

-  0.01  61eC€)*X  4  1 . 1 4  €  e i 

RETURN 

CONTINUE 

V  =  7LCC  (T) 

X  =  13. 61551 05 5756/Y  -  1,0 

511  =  (((((( (4 .9£CF-4*X  -  2.215  1F-3  )  *  X  4  7.M/3E- 

-  2.«4  17E-3)*X  4  3  .34  62E-7) «>  -  0.01825?5)*X 
4  3.61E1E-3  )*X  4  1.1  74  7922  )  *T*Y 

RETURN 

ENC 


3)  *X 


50 


4  N 


FROG  F £K  SFCTFAl  < C L  IF 17 , TAF E l , T A PF 1  ,  T flPp2 ,T A F E5 ) 

CIMENSI  (N  LSC(2e)»LCC<2S)»RCF(25),PV(25)»LPC<25(»LPL(?5),ETR(2! 
CTMFNSI  (N  CCN( 2C1,10> , 7ENP  f 2 tl> ,FROC(?01) 

CIHFNct  CN  AVTEIMF  (t  C,2E  >,  EVTEFF(40,2S) 

CIMENSI CN  CCLCNT  <l C  ,2  5  )  .HOFCNT (4C,2  5) 

rruckcr  r»  c«n  I  c  r  \  c«rc#  cm  n  *pr  i  cr>» 


1,  NBANC 

I )  LPC (I) »  LBl  (I)  ,  STF  (I ) 


^:f^^E^si  c^  racai  gci  ♦  RArc(  sc>,racc('5g 

CCNMCN  /ISFCCM/  ISF 
READ  (1,12)  TE  N  F 
CC  50  C  NCASES  =1,2 
FEWINC  5 

READ  <1  ,6)  NSPTGT 
CC  4 0  C  N  S F  =  1,  NSFTCT 
REAP  (1,6)  TEP 
FFflt  (1,41)  NLEVEL  ,NEANC  , R ^ 

CC  21  C  1=1,  NLFVEL 
READ  (1,411  ISC(  I) ,LCC  (I) , Ff) (I) ,OV (I) 
CCNT1NU  E 

C£eI£C(1  ,41?’lp£(I) ,LBL  (I) , STF (I) 
CONTINUE 

R  p  AO  (1.6)  NALT.NRLEV 
CC  221  N*  =  1.  NALT 

FEAC  (1,42)  (CCN(NZ,Nl),  NL  =  1,  NLEVE 
CONTINUE 

CONFUTE  N E AN  TENFERAURES  ANC  CCLUFN 
K  =  N°L  E V 

CC  125  NL  =  1,  NLEVEL 
CC  112  I  =  1,  NALT 
FROC(I)  =  CCN(I,NL) *TFNF(I) 

CONTINUE 
TP  =  0.  ( 

TC  =  0.  ( 

N  A  =  NALT 
CC  1 2 C  I  =  1,  K 

TF  =  TP  ♦  SINP(PRCC(NA  -  2  )  ,  1 . 0  p  ♦  5  ,  ?> 


,  NLEVEL) 
CCLUFN  C CUNT S 


C  1  2  G  I  =1  K 

TF  =  TP  4  sfNP(PRCC(NA  -  2),  1. 0^*5, 2) 

♦  SIMF(FRCC(NA  •  3>,1.0E4S,3)  ♦  SIFP(PRCC(N£  -  5 
TO  =  tC  ,  EINP(CCN(NA  -  2,Nl>  ,1.0E*E,?) 

♦  STNF  (CCN(NA  -  3,Nl),l.GE«E»?)  ♦  STNP(CON(NA  - 
CCL CNT ( K  -  I  ♦  1  ,NL)  =  TC 

TF/TC 


CCLCNf ( K  -  I  ♦  1  ,NL)  =  TC 
AVTEFPd  -  1  ♦  1,NL)  =  TF/TC 
NA  =  NA  -  5 
F CRCNT (  I , M )  =  C.C 
CONTINUE 

rease  =  e.43iP48 

CW2?  5.c  K 

TCH  =0.0 
NA  =  NALT 
RP2  =  RFASE’RBASE 

KK  r  K  -  I  ♦  1 
CC  125  J  =  1,  KK 
RTCP  =  RTCT  -  1 • CF ♦ 5 
STH  =  S  CF  T (PTC  P  *  R  T C f  -  F  E 2 ) 
HE  =  SVF  -  STH 
SVH  =  STF 

FT  OP  =  RTCF  -  1  •  CF  4  5 

STH  =  S CRT (PTOP* F TCP  -  R E 2 ) 

F4  =  SV F  -  STH 

SVH  =  STF 

Rtcp  =  FTCF  -  1 . CP*5 

STH  =  SCRT  (RTOP*RTCF  -  RE?) 

H  T  =  SV  F  -  STH 
SVH  =  STF 


5). 1.CE45,?) 

-  5 ,NL) » 1 . OF*  5 


RTCP  =  FTCF  -  _ ,  ^ 

STH  =  SCfiT (PT0P"RTcR  -  RE?) 

F  2  =  SV  F  -  STH 
SVH  =  STF 

RTOP  =  FTCF  -  1  •  CE  ♦  5 

STH  =  SCRT  (AES (RTCF*R TCF  -  RE?)) 

FI  =  SV  F  -  STH 
SVH  =  STF 

TFH  =  TFH  »  Cl  NT  (FRCCTNA  -  2>,H4,H5,2) 

♦  CINT  (FRCC (N  A  -  ’)tF!,F4.3>  ♦  QTN  i  ( FRCC  (N A  -  5),H1,F2,3> 

TCH  =  T CF  4  CINT(CCN(N£  -  2 , NL  )  ,H4  ,  F5  ,  Z> 

4  0  INT  (CCN  (  NA  -  3,NL  )  ,F3,H<.  ,3)  4  QINT(CCMN£  -  5  ,NL  )  ,H  1,  F  2  ,3  ) 

NA  =  NA  -  5 
CONTINUE 

FORCNT(I,ND  =  TCF  51 

EVTENP(I.NL)  =  TFF/TCF 


1  •  C  E  4  p 
iP»RTC^ 


*  \ 


TT 


1?» 

12S 


;3n 


*,  c  c 

E  cn 

6 

P 

1 1* 

IS 

U 1 
4? 
4? 


1CP 

no 

12  0 


ICC 

110 

120 


FPASE  =  REASE  ♦  E  .  C  £♦  S 

CONTINUE 

CCNTINU  E 

CC  4 0  C  NE  =  1,  N5ANC 
CC  2TC  I  =  1,  NPIEV 

ccn?inu2a?>  HT*Pt,:C(I)  <i>,radb  a> 

NRiTE  (4)  IL<C(I) , 1=1,251 
NRITE  (4)  <  CNCRCNTC I,J),J=1,25)  ,1=1,40) 
KRITEUM<EVTEmP<I,J),J=1,2E),I=1,mO> 
CONTINUE 
CCNTINU  E 
STCF 

FORMAT  ( El  6 ) 

F  C  R  M  A  T  C1FEE11.T) 

FORMAT  €1X1 
FCOMAT  (1CEE12.4) 

FORMAT  (INI) 

FORMAT  (216,5512.4) 

FORMAT  (8F1C.4) 

FORMAT  (F£.1,TE17.4) 

END 


FUNCTION  SIMP  (A  ,N  ,  K ) 

CIMFNSICN  A  (  7) 

IF  (  K  -  2)  130,  11C,12C 

continue 

SIMP  =  ( E  •  0  *  A  ( 2)  ♦5,C«ttll  -  A(7))*H/12,0 

RETURN 

CCNTINU  E 

SIMF  =  (  E  .  C*  A ( 2  I  4  5. C* A (31  -  A(1))*H/12.0 

RETURN 

m^u\k  .C*A(?)  4  All)  ♦  A(3))*H/3.0 

RETURN 

ENG 


FUNCTION  CINT  (A,N1,N2,K) 
CIMFNSICN  A ( 3 ) 

NS  =  HI  4  N  2 

IF  (K  -  2)  ICO  ,  1 1C, 12C 

CONTINUE 

CINT  =  (  (  A  (2) *H  S  -  All ) *H2 

♦  (A  (1  )  ♦  A  12)  I *C.  El *N1 
RETURN 

CCNT  INU  E 

GINT  =  ((( <2  >*HS  -  A  ( 1  )*H? 

♦  (  A  f  2  )  4  At?) )*C,5)*N2 
RETURN 

CONTINUE 

CINT  =  (A  (2)*HS  -  A  (1  )  *  N  2  - 
/  (E  .0*  NS*N 1*H  2)  4  (A ( 1 ) *H 1 
RETURN 
END 


A (T)*H1) »H1/ (6 ,C*H?*NS  ) 


A (T)*H1) *H2/ (6 . C*Hl *HS  ) 


A ( 3 ) *H  1  )* (Hl**3  4  H?»*3) 
4  A  (?)  *  NS  4  A  (3)  *N2>  *  P.  E 


f 


1 


52 


PPOGRAM  SFCTRA2  (CUTPU T»TAFEA»TAPE1,TAPE2,TAFE5) 
DIMENSION  LEC(?5)  ,L*C(25) 

ClMENSlCN  CCN(10I  ,T£MPJZ01) 

CIMFNSTCN  MCRCNT  <4C, <51  ,FVTEKF(4(»25) 

EiKK^sKSAJ  5&’R'CfM  50,’Rflcc<  501 

CCMMCN  NVNLM(25  0G>  ,  ENL  S  (  25  CO  >  ,S  TPL  (  £5  001  ,  HR  AC  (25  00) 
CCMMCN  SFEClieOll 
CC  6  0  C  J  =  1  *  14 
REWINC  1 
FEWINC  A 

READ  (1,12)  TEMP 
CC  500  NCASES  =1,2 
REWIND  5 

CC  10E  1*1,  1801 
SPEC(I)  =  C.O 
CCNT1NUE 

READ  (1  ,6 )  NSPTCT 
CC  4 0 C  N  5 F  =  1,  NSFTCT 
READ  (1  ,6)  T  SD 
FEAC  ( 1  ,41)  NLEVEL.NEANC.EM 

REa8  (A)  (  F  Nt  1  f  ,JKK=1.25)  ,1  =  1, AO) 

REAO(A)f  (EVTEMP ( I  ,  jKK)  ,jKK=l,<5)  ,1=1  ,4  01 
CC  21 C  1=1,  NLEVEL 
READ  (1,11)  LSCC  ,LCC  ,  RCN ,  EV 
CONTINUE 

CC  22 C  1  =  1,  N  E  A  N  C 
READ  (1,41)  L EC ( I )  ,  LEU  ,  S TR 
CONTINUE 

READ  (1,6)  NALT.NRLEV 
CC  221  N  2  =  1,  NALT 
READ  (1,4  2)  (CCN ( Nl)  ,NL  =  1, NLEVEL) 

CCNTINU  E 

IF  (ISP  .NE.  3)  GC  TC  225 
CNTHR  =  2.0*HOPCM(j,1) 

RHkVFR81t{E^6L»f?,T*PE5 

CALL  WATSTR  (RACC(j),RN,FVTENF(j,1)  ,CNTHP,1, LINES) 
CALL  FILTER  (SPE C  ,HR AC  ,KVMM  .LINES) 


^nUcun  E  =  1,NEANC 
CC  23C  1=1  ,NRLEV 


RE  AC ( 1 ,4  2  )  FT,PACC(I),RACA(I)  »R ACE ( I  ) 

<20  CONTINUE 

n5=kcccnec/igc  ,10C) 

NF=NCC  CN EC  ,100) 

IF  ((N<3  .GT.  10)  « AN  C  •  (ISF  .EC.  2))  GO  TC  40C 

CALL  ROTATE  (5, LINES) 

CALL  WATSTR  <R A CC ( J ) , F M , 8V  IF M F ( J ,NF ) , ? . 0 *HOP CNT < J ,N F ) , 2 ,LI NES ) 
CALL  FILTER. (SPEC, HRfC,WVNLN, LINES) 

A  0  0  CCNTINU  E 

WRITE  (  *,16) 

NR IT  E  (2,8)  SPEC 
500  CONTINUE 


ECO  CONTINUE 
STOP 

6  FCPMCT 
8  FORMAT 
14  FORMAT 
12  FORMAT 
IE  FORMAT 
A 1  FORMAT 
A?  FORMAT 
42  FCFMAT 
END 


(€16) 

(1FEF11.3  ) 
(IX) 

(1F6E12.4) 

(1H) 

( 21 E , 5F 1 2  .  A  ) 

(EEic.t) 

(FE  .1  ,-»El  7.4) 


SUBROUTINE  N AtST f  (RAC, RM,T,CFN,N, LINES) 

CCMMCN  /ISFCCM/  ISP 

CCMMCN  NVNLM(2500),ENLS(2FOC),STR(2500),HRAC(2500) 
TCFAC  =  <256 , C/T  )*M.  ' 

IF  (ISF  .EC.  2)  TCFAC  =  ZEE.C/T 

IF  (Jpf  .EC.  h)  TCFAC  =  2SE.C/T 

GO  TC  (  1  C  C , 2  CO , 3CC) ,  N 
C CNT  ]Nlj  E 

(S8ME  : 

:c  lie  I  -  1  •  I  INFS 


>  V  \ 


110 

260 


MHCLAM  =  1.0E44/NVNLM(I) 

TAU  =  TALE AC*STR  (I)*EXP(-1.43e79*ENLSfI)/T) 

1  RS^i0=-?!J£J?4Uf;s3^'W2ir'T'” 

FFAC(I>  =  RATE  AC* FCTE A*SFUN  (TAU) 

*FXF(-  1.43879*  (  1 . C E ♦ 4 /HHCL A M  4  ENIS<I)>/T) 
- NU  E 


CCNTIM 

RETURN 


:io 


:2o 

:G0 


CONTINUE 

TAUFAC  =  4.8815E-18*SCET  (RM/T)*CEN 
TEFAC  =  <T  -  296. C)/<T*29C. 0*1.47979 
SUM  =  0  .C 

CC  21C  1=1.  LIKES 

FSTR  =  TGFAC’STK II)*EXP (TEFAC*ENLS (I)) 

FHOLAM  =  1 • GE*4 / XVNLH ( I> 

TAU  =  T ALFAC*HHCLAM*HSTR 

HRAC(I)  =  TAL*SFLK (TAl ) *  E  XF  < *14387. S/(HHCLAM*T) ) 
/HHCL  A  M  *  *  4 
SUM  =  SIM  4  HRAD(I> 

Continue 

RAOFAC  =  FAC/SUM 
CC  22 C  I  =  1,  LIKES 
FFAr(I)  =  HRAn ( I  )*R  ACF  AC 
CCNTINU  E 
FFTL'RN 
CONTINUE 
RETURN 
EKC 


FUNCTION  SPUN  (TAL) 

?F~(TP! Jl!Ll.C)  GO  TC  110 
X  =  (T  -  2.0*0. 5 

SPUN  =  (  <  f  ( (  0.  00  106  1*X  -  0.0C4288>*>  4  G.  013686)*X  -  C. 02491  ) 

1  *X  4  0 .114932  ) *>  -  0.2E7057)*X  4  0  .556489 
RETURN 

110  CONTINUE 

IF  (I  .GT.  1C.0)  GO  TC  12C 

X  *  (T  -  7.0/3. C  „  , 

S FUN  =  (((((C.  00C3t7*X  -  0.0011E7)*X  ♦  0.003H8)*X  -  C.C05545) 

1  *X  4  0 • C28248) *  X  -  O.C81972)*X  ♦  0  .247c00 
RETURN 

120  CONTINUE 

Y  =  ALCG(l) 

X  =  1.0 /Y 

S  FUN  =  (  (  (  ((  C.  CE  787E*X  -  0.2€PS58>*>  4  C.  2645€8>*X 

1  -  0.2691E7)*X  4  0 . 32*  368  ) *X  4  1  .  1 2 83  79 ) *S0 RT (Y ) / T 

RETURN 
END 


SUBROUTINE  FILTER  (CPEC  .RAC.NVNUK.N  ) 

iKi£,""a 

COMMON  / I S  RC  CM /  ISF 
CC  1 5 C  I  =  1,  N 


XV  N 


mill' 


120 

150 


RT  = 

RADT 

=  (PT  -  4 0*> 
JF  =  (RT  -  395 
IE  (JS  .IT.  1) 
IF  ( jF  .OT 
IF  („S  .GT 
CC  1 2 C  J  = 

R  N  U  M  =  (J 
SPEC  ( j)  = 
CONTINUE 
CONTINUE 
R  R  T  U  R  N 
END 


0 )*C.5  4  2  .  C 
0  >*0.5  4  1,0 

JS  =  1 

1801)  jF  =  18C1 
.  JF)  GO  TC  1*0 
JS.JF 

•  1)  *2  .0  4  L  CC.O 
SPEC(J)  ♦  ( 5  .  C  -  AES  (RNL'M 


PT)  )  *  R  A  CT*  0.0  4 


SUSP CUT  INF  ROTATE  (NlN, LINES) 

CCMMCN  n'FCCM/  ISF 

COMMON  NVNLM  (2  500, ENIS(  25  CO), STP(25C0),HFAC(  25  00) 
P  E  AC  (NIN)  LINES 

IF  ( EOF ( NI N  )  )  121, 19C 

190  READ  (NIN) (KVNUM? I) ,STR  (I) ,ENLS( T) ,1=1. LINES) 

IF  (ECC( NIN  )) 121 , 115 
119  CONTINUE 
121  RETURN 
ENC 


c 


FivOGhfl‘1  FLCTKt  I  N  FUT,CUTPL,T  =  /PCP0,TAFE1  =  INPUT,TARE2=CLTPLT,  TAPED 
LAST  CHANGE  PER  IE,  197° 


1 

5  0  C 

1C00 
1C 
55C 
1  9 

P63 

?  ccon 

^OCdO 


-.0000 

80 


«cc 

9  c  n 


DI  NFN  SION  *<ie) 

DATA  COL  /•’?"/ 

FFWTNC  3 

R£AD(3,5CC>  TY°E 
FCrMAT(A  1.0  7,218, 7^18. 11) 
TF  (  T  Yc  E.  EC  .COL)  GC  TO  1C 
REWIND  3 


WRITE!?,  100C) 

FCFNAH*  E N T E K  CELT, START, STCF,ANP, OFFSET, STAPTO*) 
REACH  1,* )  CELT, WTNC, XENC«ANF  ,CFF, START 

READ(3,9  00)  T  YD£  ,  TITLE  , NUN ,  Ik E V , AF , D EL T , ?T APT 
WRITE (2, 950  ) 

ENT-p  START  .STCF.ANF, OFFSET*) 

°  £  00 ( 1 ♦ ^  I  WIND, XENC, APF  ,CFP 

P.EKINC  2 

xnELT=PSC./tX  END- WIND) 

XDFLf =tplt* XTELt 

wf:te(2,eeo)  title 

FCMATI/  ,15X,  AT  ) 

PORNaH^X^IN^^VAL^E^WE^M 
W^I TE(2,3000 C)  CFF,ANF 
FCRMAT!3>,*ANPlITLC£  RANGE!*, F?.0,*  TO  *,F7.0,*)*) 

RFKlNC  2 
CAlL  FLOTS 
CALL  PRANE 


,  F  8 . 3  ,*  AND  *  ,FP.3> 


IY2=8 

1X2=10 

X?=10 

XN=  STAR! 


CO  ICO  1=1 , l C0000 
KEAC(3,LOOOO)(X(J),J=1,1E) 
FORNAT (16T5) 

Ic ( E  CP ( 3 ) )  9CC,80,?0 

X  =  X  N 

XN=  >M16*DELT 
IF( XN.LT.WINCJGCTC  100 
CO  EC  J=  1 , 1  6 
XJ=X4FLOAT( J-l) *CELT 
IP(  XJ.LT »WINC)GCTC  90 


IX= IX? 

IY= IY? 

X’ =X24X0FLT 
IX2  =  X  2 

IP (IX2.GT.9C0IGCTC  800 
:  Y  2=5  O’J  •  *(<(J)-CFP>  /  (ANP-CFF)  4®, 5 
IF ( IY2 , L£  •  1 )  I Y  2  =  1 

CALL  LINE!  IX  ,IY,IY2,IY2) 

mn  mi 


REWIND  3 
CALL  XNTT 
PP  AD !  1,* ) JUNK 
REWIND  2 
GOT  C  1 
5TOF  7 
End 
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ooo  !-*■  m  ooo  no 


r  THIS  PTUTINE  INITIALIZES  thf  PLOT  EUFFEP 

C  ANC  THE  INCEXFS  INTO  THE  PLOT  EUFFEP • 

r, 

SLDPCUTHF  °L  CTS 

C3MMCN/IXY/I T,irX,ICLY,ICHY,TCLX, IOHX 
COMHCN/DUF  FFR/LB, LK,IB(5C0» 


L  E  =  1 
IT  =  0 
L  t*  =  3  f 

IF(l»=C01TOOCCCOCOOOOflOOOCl3 

RETURN 

eND 


THT  S  SUBROUTINE  PACKS  A  12-RIT 
PLOT  IMSTPUCT ICN  FCR  THE  TEKTRONIX  INTO 
T Hf  FUFFE?  (IB).  IT  IE  FCPMA7TEG  TO  THE 
PECULIAR  BIT  S  T  Fr  AM  F EGUI FEMEN'TS  FCR  TcMFO • 

SUFPO  ITTNE  SENE  (I) 

CCNMON/DUFFER/lBiLKtlRtiOO) 

IF  (LW  .GE.OGOTC  2 
LW=u* 

L a  =LB  -*1 

IB (LB  )=CR ( ANO (  IE (LB)  ,SHIFT  (^ASK(4  f) »LW)  ),SHIFT(I,LW) ) 
LW=  LN-19 

FOPMA  T  (  3X,I4,  3X,T4) 

PrT LRN 
END 


1 

1 


CCMM  CN  /  IY  Y/  ITt^CX*ICLY»ICH>*IOLX  »I  OHX 
CCMMCN/EUFFER/  LE.LW  ,IP(c00) 

BStt  IKEliKTj 

CALL  S  ENC  (TOHX) 

CALL  S  ENC(IOLX) 

DC  1  K=l.° 

IF  (LW.EC.  CTGOTO  2 
CALLSENC(  C) 

BUFFER  CUT  (?,l)  )  (  IE  (  1 )  ,  18  (  LB)  > 


CALL  FLCTS 

RETURN 

END 


V/ 


c 


SLDRCUTINE  pacf 


CALL 


CALL  SENC(403ZE) 
CALL  SENCUO  1  4B) 
CALL  SENG (4  00  7F) 
SEN  C  ( C) 

CALL  xv.it 
P  ETU  P  N 


END 


SUFFOITINE  COPV 


CALL  SENC<40  3^E) 

Ha bh  3»?c*0?7n 

RETUR  K 
ENC 


SL:JP  CLTINE  FRAME 
CALL  L  IRE  (10,  50  9,900 ,509) 


CAL1- 

CALL 


;all 

CALL 
CALL 
CALL 
CALL 

call 

CALL 
CALL 
CALL 
CALL 
CALL 
PPT  URN 
END 


LINE( 9,  rfa  8  *  9 , 9  ) 

LIKE( 898 , 7* 8,898,9) 
LINE  <  1  0,  7fc9,90C,  76  8) 
LlNEIlO,  767,90C«  767) 
LIKE(10  ,10,900,10) 
LlNEIlO, 9,900,9) 

LT  KP (1  0  ,  76  8,1  0,9) 
LlRE(li.7f 8,11,9) 
LlNE<9  00,76«,900,9> 
LIKE (999, 76 8, 899,9) 


SUERCLTINE  L IKE C J  ,  I , J J , I T ) 
CCMMCK/EL'FFER/LP 
CALL  SEN  C(40T58) 


K= 09 ( A  KCC SHIFT (1,-51 , -MASK ( 95 >  > , 40  10") 
CALL  9  ENC (K) 

K=OR(AKC(I,-FASK(55)),A140E» 

CALL  SENC(K) 

K  =  OR(AND (SHIP  T( J,-5)  ,-KASK(55)>,40l0B) 
CALL  S  E  N  C  ( K  ) 

K=OF(AKC(j,-VASK(55)),4inOP) 

CALL  S  ENC  (K) 

K=OR<ANC  (SHIFT  (11,-5)  ,-rMASK  (95))  ,4C40P) 
CALL  SENC(K) 

K=00< ARC (11,-MASK (55) )  ,4140E) 

CALL  SENC(K) 

K  =  OR(AKC(SHIFT ( JJ,-5)  ,-MASK  (55) ) ,4C40B) 
CALL  SEKH(K) 

k=Cp<ANC(jJ,-MASK(55)),41C0E) 

CALL  SENC(K) 

IF(lB  .GT.«,80)C  ALL  XMTT 

RETLFN 

END 
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Appendix  C 


• 


PLOTR 

EQU 

GF5H 

• 

PLTCS 

EQU 

0F7H 

• 

MOMCS 

EQU 

OEFH 

MDMQT 

EQU 

DEEM 

TXROY 

EQIJ 

Gi 

RXRDY 

EQU 

2 

8HODO 

EQU 

83H 

EOCHR 

EQU 

83H 

• 

* 

NDCHR 

EQU 

0  FFH 

LFTPN 

EQU 

43  H 

• 

ORPEN 

EQU 

8  OH 

• 

t 

LWAIT 

EQU 

50 

• 

« 

TIME  DELAY  FOR  PEN  UP/DOWN 

SWA  IT 

EQU 

255 

OVRIT 

EQU 

10H 

0  ICMP 

EQU 

1QQ1Q0Q1B 

• 

* 

NXORG 

EQU 

0600 

• 

1 

•  ** 

:  »* 

ORG 

38  0  0  H 

IN 

NDMDT 

• 

• 

CLEARS  MOOEM 

M  VI 

A,  *R ' 

• 

CYBER  RNH  COMMAND 

CALL 

XMIT 

MVI 

A,  *N  * 

CALL 

XMIT 

MVI 

A,  *H  * 

CALL 

XMIT 

CALL 

SNDCR 

• 

XMITS  • RETURN  * 

CALL 

HEXIN 

• 

* 

IGNORE  CR,LF:WAIT  FOR  PROMPT 

MVI 

A,  BMOOC 

• 

I NIT  PLOTTER  PORT 

SET  55 

EQU 

S 

OUT 

PLTCS 

;♦* 

L  X I 

H,  BUFER 

• 

* 

RESET  PNTER 

SHLO 

PM  TER 

:♦* 

L  X I 

H,  0 

• 

* 

CLFAR  CURPX,  CURRY 

SHLO 

CURRX 

SHLD 

CURPY 

:** 

MNITR 

EQU 

t 

CALL 

NXCHR 

• 

• 

GET  COMMAND 

ADD 

A 

• 

• 

OFFSET  =  2*COMMAND 

MOV 

C.  A 

MVI 

B,  0 

!*• 

L  X I 

H,  COMNO 

• 

• 

GFT  ADDRESS  OF 

OAO 

3 

• 

• 

command  routine 

MOV 

- »  M 

• 

• 

INTO  DE,  THEN 

I N  X 

H 

MOV 

XCHG 

Of  M 

• 

INTO  HL,  THEN 

PCHL 

• 

t 

INTO  PC 

:•* 

COMNO 

EQU 

5 

OW 

pnlft 

• 

t 

MOVE  W/  PEN  UP 

QW 

OW 

pndrp 

XGRAF 

• 

• 

• 

« 

MOVE  W /  PEN  DOWN 
X-AXIS  GRAPH  ROUTINE 

DH 

YGRAF 

• 

• 

Y-AXIS  GRAPH  ROUTINE 

OH 

TXPLT 

• 

• 

PLOT  CHARACTERS 

DW 

RETCR 

• 

♦ 

RETURN  CURR.  X.Y 

OW 

SETCF 

• 

* 

RESET  CURR.  X,Y 

cw 

ENDIT 

• 

1 

TERMINATES  PLOTTING 

:«* 

• 

* 

COMMAND 

ROUTINES 
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PNLFT  EQU  $ 

CALL  LIFT 
JMP  PNMOV 

PNORP  EQU  S 

CALL  DROP 

PNMOV  EQU  * 

CALL  REGIN 
SHLD  ENDX 
CALL  REGIN 
SHLO  ENOY 

CALL  PLOTS 

J MP  MNITR 

•** 

!** 

XGRAF  EQU  $ 

CALL  NXCHR 
MOV  C,  A 
LHLO  CURRX 
SHLO  ENDX 
CALL  REGIN 
SHLO  ENOY 
CALL  LIFT 
CALL  PLOTS 
LOOPX  EQU  $ 

CALL  REGIN 
MOV  A,  H 
CPI  NOCHR 
Jl  MNITF 

SHLO  ENOY 

LHLO  CURRX 
DAO  B 
SHLO  ENDX 

CALL  DROP 
CALL  PLOTS 

JMP  LOOPX 


YGRAF  EQU  S 

CALL  NXCHR 
MOV  C,  A 
LHLO  CURRX 
SHLO  ENDX 
CALL  REGIN 
SHLO  ENOY 


CALL  LIFT 
CALL  PLOTS 
LOOPY  EQU  % 

CALL  REGIN 
MOV  A,  H 
CPI  NOCHR 
JZ  MNITR 

SHLO  ENOY 

LHLO  CURRX 
OAO  3 
SHLO  ENOX 

CALL  OROP 
CALL  PLOTS 

JMP  LOOPY 


TXPLT  EQU  % 

CALL  NXCHR 


:*♦ 

:♦* 

•  £ 

:  LIFTS  PEN  IF  NEC 

:  * 

*  # 

1  DROPS  PEN  IF  NEC 

!  # 

:  # 


5  READS  ENDPTS  ANO 

:♦* 

;  MOVES  PEN  THERE 

s  i 


PLOT  X-AXIS  GRAPH 
GET  X-INCR. 


B  IS  CLEAR 


MOVE  PEN  TO  1ST  POINT 


;  PLOT  STREAM 
J  GET  NEXT  Y  POINT 
%  CHECK  FOR  ENO  OF  STREAM 


:** 

?  UPDATE  X  COORD 


;  INSURE  PEN  DOWN 

;+• 

;*♦ 

5  PLOT  Y-AXIS  GRAPH 
J  GET  Y-INCR. 

?  B  IS  CLEAR 
:  MOVE  PEN  TO  1ST  POINT 


!  PLOT  STREAM 
?  GET  NEXT  Y  POINT 
!  CHECK  FOR  ENO  0F  STREAM 


!*♦ 

!** 

;  UPDATE  X  COORD 


;  INSURE  PEN  DOWN 

1  ** 

;  GET  SCALE 


.  .«■  v\ 


CLOOP 

A  DOFF 

GTSEG 

SLOOP 

RETCR 

SFTCR 

ENOIT 

REGIN 


MOV 

E,  A 

• 

4 

AND  SAVE  IN  E 

;•* 

CALL 

NXCHR 

• 

4 

GE.  '•IRECTION 

STA 

OIRN 

!»• 

•** 

EQU 

% 

CALL 

NXCHR 

• 

♦ 

GET  NEXT  CHAR.  CODE 

CPI 

NDCHR 

• 

4 

RETURN  TO  MONITOR 

JZ 

MNITR 

• 

4 

AT  ENC  OF  STRING 

LHLD 

CJRRX 

• 

♦ 

RESET  ORIGINS 

SHLD 

ORIGX 

• 

9 

* 

LHLD 

CURRY 

• 

• 

# 

SHLD 

ORIGY 

• 

4 

# 

LXI 

H.  CHTAB 

• 

4 

DATA  TABLE  FOR  CHAR  GEN 

ADD 

A 

• 

4 

OFFSET=CODE*2 

JNC 

A30PF 

!»♦ 

INK 

H 

EQU 

l 

AOO 

L 

• 

• 

A  CONTAINS  OFFSET 

MOV 

L.  A 

• 

9 

♦ 

JNC 

GTSEG 

4 

4 

♦  HL  POINTS  INTO  CHTAB 

I  NR 

H 

4 

4 

# 

EQU 

$ 

MOV 

C*  M 

• 

♦ 

BC  POINTS  TO  SESS 

I  NX 

H 

MOV 

R  *  M 

EQU 

« 

• 

• 

« 

LOAX 

3 

• 

9 

GET  NEXT  SEGMENT  CODE 

MOV 

0,  A 

• 

4 

SAVE  SEG.  CODE  IN  0 

CALL 

CHGEN 

• 

4 

PLOT  SEGMENT 

I  NX 

8 

• 

4 

INDEX  OF  NEXT  SEGMENT 

MOV 

A,  C 

• 

♦ 

# 

ANI 

EOCHR 

• 

4 

♦CHECK  IF  lAST  SEGMENT 

JZ 

SLOOP 

• 

4 

CONT  IF  NOT  LAST  SEG 

MVI 

At  NXORG 

• 

• 

♦PEN  UP,  (u,6»  *  SCALE 

MOV 

D »  A 

CALL 

CHGEN 

• 

4 

INTERLETTER  SPACE 

JMP 

CLOOP 

4 

4 

F CP  NEXT  CHAR 

EQU 

S 

4 

LHLD 

CURRX  ; 

SEND  CURRX.CURRY 

CALL 

PGOUT 

LHLO 

CURRY 

CALL 

RGOUT 

;** 

CALL 

SNOCE  T 

SEND  'RET  * 

CALL 

HEX  IN  ? 

WAIT  POP  PROMPT 

JMP 

MNITF 

t*» 

EQU 

t 

CALL 

REGIN 

SHLD 

CURRX 

CALL 

REGIN 

S  HI  0 

CURPY 

JMP 

MNITP 

;*• 

:** 

EQU 

It 

RST 

1  ? 

RETURN  TO  MONITOR 

:** 

•  ** 

• 

4 

I/O  UTILITIES 

:** 

EQU 

E  ; 

4  HEX  DIGITS  INTO  HL 

CALL 

NXCHR 

MOV 

H,  A 

CALL 

NXCHF 

MOV 

L,  A 

RET 

:»* 


NXCHR 


LOONE 


LOADR 


LODNG 


HEX  IN 


EQU  t 
PUSH  B 
PUSH  0 
PUSH  H 
LHLO  PNTER 
LX  I  D,  BUFEk 
CALL  CMP  AF 
CZ  LOADR 

MOV  9,  M 
I  NX  H 

PUSH  H 

LXI  Ot  SUPER 
LHLD  BUFLN 
DAO  0 
XCHG 
POP  H 

CALL  CMPAR 
JN  7  LOONE 

LXI  H,  BUFER 
EQU  S 
SHLD  PNTEF 
MOV  A,  B 
POP  H 
POP  D 
POP  9 
RET 


1 

:** 

;  GETS  HEX  CONSTANT  from  BUFFER 

:  CHECK  IF  BUFFER  EMPTY 

5  FILL  IT  IF  NECESSARY 

:  ** 

\  GET  CHARACTER 

:•* 


:** 


ELSE 


RESET  PNTER 


EQU 

* 

CALL 

LIPT 

PUSH 

D 

PUSH 

H 

LXI 

H,  Bl 

CALL 

SHAKE 

CALL 

HEXIN 

MOV 

D,  A 

CALL 

HEXIN 

MOV 

E  *  A 

XCHG 

• 

• 

SHLD 

BUFLN 

XCHG 

EQU 

S 

CALL 

HEXIN 

MOV 

M»  A 

I  NX 

H 

OCX 

D 

XR  A 

A 

ORA 

E 

JN  7. 

LODNG 

ORA 

0 

JN  Z 

LQONG 

CALL 

HEXIN 

POP 

H 

POP 

D 

RET 

EQU 

S 

PUSH 

B 

CALL 

A5CHX 

RLC 

RLC 

RLC 

RLC 

MOV 

CALL 

aIchx 

M* 

:»» 

;  PREVENTS  INK  SPOTS 

:** 

?  GET  WORC  COUNT 

SAVE  BUFFER  LENGTH 


?  GET  NEXT  WORO 


\  CHECK  IF  LAST  WORD 


\  REPEAT  UNTIL  DONE 
:  WAIT  FOR  PROMPT 

!»• 


;♦* 

!** 

;** 

?•* 

;  2  HEX  DIGITS  INTO  H 
;  ASCII  TO  HEX  CONVERSION 


SAVE  HIGH  CHARACTER 
GET  LOW  CHARACTER 


6l 


ASCHX 

EQU 

s 

* 

* 

i** 

ASCII  TO  HEX  IN  A 

CALL 

RECEV 

• 

A 

GET  ASCII  CHAR 

CPI 

JC 

•A* 

ARNOA 

• 

« 

JUMP  IF  NUMBER 

AOI 

9 

* 

♦ 

:♦* 

ADJUST  A-F 

ARNOA 

EQU 

ANI 

$ 

ofh 

• 

A 

CLEAR  HIGH  BITS 

RET 

!** 

:♦* 

RECEV  EQU 

* 

MDMCS 

ANI 

RXROY 

JZ 

RECEV 

IN 

MDMOT 

ORA 

A 

JZ 

RECEV 

CPI 

JAM 

JZ 

RECEV 

CPI 

QOH 

JZ 

RECEV 

RET 

RGOUT  EQU 

$ 

CALL 

HXOUT 

MOV 

H.  L 

CALL 

HXOUT 

RET 

HXOUT  EQU 

* 

MOV 

A.  H 

RRC 

RRC 

RRC 

RRC 

CALL 

HXASC 

MOV 

At  H 

CALL 

HXASC 

RET 

HXASC 

EQU 

1 

ANI 

OFH 

CPI 

0  AH 

JC 

ARNOC 

AOI 

7 

ARNOC 

EQU 

S 

AOI 

•0  • 

CALL 

RET 

XMIT 

XMIT 

EQU 

$ 

PUSH 

9 

MOV 

Ct  A 

XLOOP 

EQU 

* 

IN 

MOMCS 

ANI 

TXROY 

JZ 

XLOOP 

TEST  MO DEM  STATS 

WAIT  FOR  CHAR 

;** 

REAC  CHAR 

IGNORE  NULL  CHARS 

•  •  •  L  F  v 

•  •  •  C  R 


:•* 

:  (HL>  TO  4  ASCII  CHAR 


:** 

;** 

!♦* 

:  (HI  TO  2  ASCII  CHARS 


:  LOW  4  BITS 


:♦* 

1  (A)  TO  1  ASCII  CHAR 
:  SEND  LOW  4  BITS 

*  CHECK  FOR  AND 

?♦* 

;  ACJUST  A  -  F 
t  FOR  ASCII  "PREFIX" 


!♦* 


MOV 

OUT 

mShot 

CALL 

RECEV 

POP 

B 

RET 

SHAKE 

EQU 

% 

PUSH 

H 

M  VI 

L.  CAH 

EWAIT 

EQU 

S 

CALL 

LONG 

OCR 

L 

JNZ 

EWAIT 

POP 

H 

M  V I 

A,  *1* 

CALL 

XMIT 

SNDCR 

EQU 

* 

IN 

MQMCS 

ANI 

TXRCY 

JZ 

SNOCR 

M  VI 

A,  QDH 

OUT 

RET 

MDMOT 

• 

« 

PEN  CONTROL 

LIFT 

EQU 

S 

LOA 

PSTAT 

ORA 

RZ 

A 

MVI 

A,  LFTPN 

CALL 

MOVE 

• 

JMP 

PWAIT 

DROP 

EQU 

5 

LQA 

PSTAT 

ORA 

RNZ 

A 

MVI 

A,  CRPrN 

CALL 

MOVE 

• 

INR 

A 

PWAIT 

EQU 

5 

STA 

PSTAT 

LONG 

EQU 

% 

PUSH 

H 

MVI 

L,  LWAIT 

LLOOP 

EQU 

? 

CALL 

SHORT 

OCR 

L 

JNZ 

LLOOP 

POP 

RET 

H 

MOVE  EQU 

S 

OUT 

PIOTP 

X  R  A 

A 

OUT 

PL  OTP 

CALL 

SHORT 

:»* 


•  «♦ 


5  WAIT  FOR  CYBER  FQ X 
»  APPROX.  1  SEC 


5  HANDSHAKE  CHAR 

:  XMTS  *  C F  * 


:  *  FFTURN  * 

:  FCHO  «.  LF  ignored 

: 

UTILITIES 

s  it 
:  it 
:  it 

:  # 

• 

i  CONTROLS  PLOTTER 
5  * 

;  # 

:  # 
j  it 
:  h 

•  «•.*****■»*♦♦ 

;  CONTROLS  PLOTTER 
;  PSTAT  GETS  1 

:  it 

:  « 

;  * 

:  it 


i  n 
:  it 
:  * 

;  ft 
:  # 

:  *■» 

;  SEND  INSRTUCTION 

;  FOLLOWED  3Y  C 
•  WAIT  1.65  M.LLISEC 


!  ** 


SHORT  EQU 
PUSH 
M  V I 

MLOOP  EQU 
OCR 
JNZ 
POP 
RET 


SWAIT 


S 

L 

WLOOP 

H 


ARITHMETIC 


:  * 

:** 

UTILITIES 


0PSU8  EQU 
MOV 
SUB 
MOV 
MOV 
SBB 
MOV 
XRA 
RET 

CMPAR  EQU 
MOV 
CMP 
RNZ 
MOV 
CMP 
RET 

SCALR  EQU 
AMI 

MULT  EQU 
PUSH 
PUSH 
M  V  I 
MVI 
MOV 
MOV 

MLOOP  EQU 
RAR 
MOV 
JNC 
0  AO 
CMC 

ARNDM  EQU 
XCHG 
DAO 
XCHG 
OCR 
JNZ 


t 

B 

0 

C«  8 


8«  A 
ARNDM 
D 


C 

MLOOP 


DOUBLE  PRECISION 
SUBTRACTION  ROUTINE 

• 

* 

* 

RESULT  IN  HL,  ACCUM  CLEARED 


OCUBLE  PRECISION  COMPARISON  ROUTIM 

* 

# HIGH  BYTES  NE  3 

* 

* 

» 


# 

s** 

8X8  MULTIPLICATION 

ROUTINE 

HL  GETS  A*E 


SHIFT  MULTIPLICAND  LEFT 

* 


MAIN-LINE  UTILITIES 
!»♦ 


CHGEN  EQU  $ 

PUSH  9 
PUSH  0 
PUSH  H 

ANI  1300 
JNZ  CMKON 

CALL  LIFT 
JMP  RELPT 

CHKDN  EQU  t 

CALL  OROP 


;** 

;  # 

:  SAVE  REGS. 


:  PEN  UP/DOWN  FIELD 
;  OECODE  PEN  UP/OOWN 

:♦* 

:  # 

:  # 


RELPT 

EQU 

S 

MOV 

CALL 

A.  0 
SCALR 

MOV 

3,  H 

• 

MOV 

C,  L 

♦ 

MOV 

A,  0 

RRC 

RRC 

RRC 

CALL 

SCALR 

• 

XCHG 

* 

LOA 

3IRN 

MOV 

H,  A 

ORA 

A 

• 

JPE 

DIRX 

PUSH 

B 

MOV 

3,  0 

MOV 

C,  E 

POP 

D 

DIRX 

EQU 

5 

ANI 

19  B 

• 

JZ 

OIRY 

XRA 

A 

SUB 

E 

MOV 

E,  A 

M  VI 

A,  0 

S  33 

D 

• 

MOV 

0.  A 

DIRY 

EQU 

$ 

MOV 

A,  H 

ANI 

31 B 

• 

JZ 

SETPT 

* 

XRA 

A 

SUB 

C 

MOV 

C.  A 

M  VI 

A,  0 

S  39 

9 

MOV 

B«  A 

SETPT 

• 

EQU 

S 

• 

LHLD 

ORIGX 

DAO 

D 

SHLO 

ENOX 

• 

* 

lhlo 

ORIGY 

DAO 

B 

• 

SHLO 

ENDY 

• 

• 

CALL 

PLOTS 

* 

POP 

H 

POP 

c 

POP 

3 

RET 

CHTAB 

EQU 

K 

DM 

CHR26 

DW 

CHR3? 

DM 

CHR34 

OW 

CHR35 

ow 

CHR36 

DW 

CHR37 

OW 

CHR38 

DW 

CHR39 

DM 

CHR40 

DW 

CHP41 

OH 

CHRL  ? 

* 

0  CONTS  SEG  CODE 
(A. AND. 7)* SCALE 
LOADS  BC  WITH 
REL-ATIVE  Y 


(E) 


;  * 

x  * 
x  # 

;  SHIFT  FOR  RELATIVE  X 
:  IGNORES  HIGH  BITS  IN  A 

X  STORE  REL  X  IN  DE 

J  COMPUTE  DIRECTION  OF  TEXT 

;  # 

?  SETS  FLAGS 

X  JUMP  IF  ON  X-AXIS  ( 0 S. 3 > 

X  SWAP  X  AND  Y  VALUES 
:  FOR  VERTICAL  LETTERS 
X  ( 112> 


A 

# 

COMP  X  VALUE 

# 

* 

# 

* 

DE  GETS  -DE 


* 

* 

# 

COMP  Y  VALUE 

* 

# 

* 

* 

BC  GETS  -BC 

# 

* 

# 


:  * 

;  # 

:  * 

:  * 

:  RESTORE  REGS. 


:*♦ 

:  # 

;  ** 

:  data  for  char  gen 

;  CHARI  I 
:  CHAR*  A 
!  CHARI  9 
S  CHAR*  0 
!  CHARI  D 
5  CHAFt  F 
:  CHARI  F 
;  CHAR*  G 
;  CHARI  H 
;  CHARI  I 
:  C HA Rt  J 


65 


ow 

CHR43 

CHAP  : 

< 

OH 

CHR.44 

CHAP  « 

L 

DW 

CHR45 

CHARS 

'I 

OW 

CHR45 

CHAPS 

N 

OW 

CWP47 

CHAk  8 

3 

ow 

CWR48 

CHAPS 

p 

OW 

CHR49 

CHAP  S 

0 

OW 

CHR5G 

CHAP  s 

R 

ow 

CHAP  S 

3 

DW 

CHR52 

CHAP  S 

r 

ow 

CHRS3 

CHAP  s 

j 

ow 

CHR54 

CHARI 

v 

ow 

CHP  55 

CHAPS 

w 

ow 

CHR56 

CHARS 

X 

DW 

CHR57 

CHAPS 

Y 

DW 

CHR58 

CHAR  S 

7 

DW 

CHP16 

CHAP  S 

j 

ow 

CHR1  7 

CHAP  I 

< 

ow 

CHRlfc 

CHAP  S 

2 

ow 

CHR19 

CHAP  S 

3 

ow 

CHF2C 

CHAPS 

4 

cw 

CHR21 

CHAR  S 

5 

ow 

CWR  1  2 

C  HA1'  S 

6 

ow 

CHP  2  3 

CHARS 

7 

0  W 

CHP24 

CHAP  5 

8 

ow 

CHR25 

CHAR  S 

3 

ow 

CHR1 1 

CHARS 

f 

ow 

CMR1  3 

CHAPS 

- 

ow 

CH?1  C 

CHAP  8 

« 

ow 

CHR15 

CHARI 

/  j  L  m  3  H 

ow 

CHP  C  8 

CHAP  S 

( 

ow 

CHR  G  9 

CHAkS 

) 

ow 

CHRC4 

CHAP  8 

i 

ow 

CHR29 

CHAP  s 

= 

ow 

CHR69 

CHAR  S 

SL  ANK 

ow 

CHR12 

CHAP  S 

* 

ow 

CHR14 

CHARS 

• 

ow 

CHRG3 

CHAP  S 

«POUNC  31 

GN 

DW 

CHR59 

C  H  A  P.  5 

[OPEN  SQ 

3KKT 

OW 

ow 

mi 

:  CHAPS 
:  chaps 

) 

PERCENT 

ow 

CHRQ2 

:  CHARS 

**  COL  •  QUOTE 

ow 

CHR63 

:  CHARS 

.UNDERSCORE 

ow 

CHRO  1 

:  CHARS 

Exclamation 

ow 

CHRC6 

;  CHARS 

s.  ampersand 

DW 

CHRG7 

;  chaps 

*  APOSTROPHIE 

DW 

CHP.31 

J  CHAPS 

7QUEST ION  MARK 

OW 

CHR23 

;  CHAPS 

< 

ow 

CHR3C 

;  chaps 

> 

OW 

CHR32 

:  CHARS 

3  COMMERCIAL  AT 

OW 

CHR6G 

5  CHARS 

\ 6 A ck  slash 

ow 

CHR62 

;  CHARS 

''CIFCUMFLEX 

ow 

CHR27 

\  CHARS 

• 

DW 

CHR64 

!  CHARS 

A  POSTEOPHIE 

ow 

CHR65 

:  CHARS 

OPEN  PAPEN 

ow 

CMR66 

;  CHARS 

VERTICAL  LINE 

ow 

CHR67 

:  CHARS 

CLOSE  PAREN 

ow 

CHR68 

;  CHARS 

:  ** 

TILOE 

DHLST 

EQU 

? 

:  CHARACTER  TABLES 

•  *» 

DHR01 

EQU 

1 

;  CHARS 

OB 

0160 

08 

1120 

08 

OHO 

08 

3100 

!»♦ 

6hR02 

EQU 

% 

:  CHARS 

08 

0140 

OB 

1260 

OB 

G  340 

03 

3460 

• 

;*• 

CHR03 

"EQU 

* 

;  CHARS 

1 

OB 

0130 

C)«*  o«» 


HR04 


HR05 


HR  0  6 


CHROr 

&HR06 


CHR09 


m 

6hrio 


DB 

1260 

OB 

0360 

08 

1230 

OB 

0  350 

00 

1150 

08 

08 

EQU 

t 

08 

0020 

OB 

1110 

DB 

1310 

08 

1420 

00 

1330 

08 

1130 

08 

1040 

OB 

1150 

08 

1350 

08 

1440 

08 

0260 

08 

3200 

EQU 

S 

D  3 

1460 

08 

1160 

81 

H?8 

08 

1250 

08 

1160 

DB 

0410 

08 

1320 

08 

1210 

08 

1300 

08 

3410 

EQU 

S 

08 

0420 

08 

1200 

08 

itoo 

08 

1010 

OB 

1120 

DB 

1350 

08 

1260 

OB 

1150 

OB 

3400 

EQU 

S 

OB 

0260 

08 

3240 

EQU 

S 

08 

0360 

08 

1240 

08 

1220 

08 

3300 

EQU 

i 

08 

0160 

08 

1240 

03 

1220 

08 

3100 

EQU 

t 

OB 

0130 

08 

1350 

03 

0250 

00 

1230 

OB 

0150 

OB 

3330 

:♦* 

:  chari  i 


:♦* 

:  CHARI  PERCENT 


S  CHARI  8, 


I  CHARI ' * 


5  CHARI  ( 


:*♦ 

;  CHARI  ) 


;  CHARI  * 


:•* 


CHAM  *■ 


CHR11 

EQU 

t 

03 

0130 

09 

1330 

03 

0240 

09 

3220 

• 

CHR12 

EQU 

S 

03 

0220 

08 

3100 

• 

3HR13 

EQU 

? 

03 

0130 

03 

3330 

CHR14 

EQU 

8 

09 

QUO 

09 

3100 

• 

CHR15 

EQU 

$ 

OB 

3460 

• 

CHR16 

EQU 

S 

OB 

01  QO 

03 

1020 

03 

1340 

03 

1160 

03 

1360 

DB 

1440 

03 

1420 

09 

1300 

03 

3100 

• 

:m«i7 

EQU 

S 

09 

0140 

08 

1260 

03 

3200 

• 

CHP18 

EQU 

S 

03 

0050 

D  9 

1160 

03 

1360 

09 

1450 

08 

1440 

OB 

1120 

03 

1300 

03 

3400 

• 

CHR19 

EQU 

% 

09 

0010 

03 

1100 

03 

1330 

03 

1410 

03 

1420 

03 

1330 

08 

1230 

OB 

0330 

03 

1440 

03 

1450 

03 

1360 

03 

1160 

03 

3050 

• 

5HR20 

EQU 

S 

03 

0300 

09 

1360 

09 

1020 

08 

34  20 

:♦* 

:  CHAR  I  « 

;♦* 

:  charj  - 

:  chari  . 

;  chap.i 

:»* 

t  CHARI  3 


:** 

:  CHARI  1 

:*♦ 

;  CHARI  2 


«** 

1  CHARI  3 


:  CHARI  4 


CHR21 

EQU 

$ 

• 

CHARI 

5 

8? 

1118 

OB 

1300 

OB 

1410 

08 

1430 

OB 

1340 

OB 

1140 

08 

1330 

08 

1060 

OB 

3460 

• 

:** 

6hR22 

EQU 

S 

• 

CHARI 

6 

OB 

0030 

OB 

1140 

OB 

1340 

08 

1430 

OB 

1410 

08 

1300 

03 

1100 

DB 

101O 

OB 

1050 

OB 

1160 

03 

1360 

DB 

3450 

:♦* 

• 

5  HR  2  3 

EQU 

S 

• 

CHARI 

7 

OB 

0060 

08 

1460 

OB 

3100 

;*♦ 

6hR24 

EQU 

3 

• 

♦ 

CHARI 

8 

OB 

3100 

OB 

1010 

OB 

1320 

DB 

1130 

08 

1040 

OB 

1050 

OB 

1160 

DB 

1360 

OB 

1450 

OB 

1440 

08 

1330 

08 

1130 

OB 

0330 

DB 

1420 

OB 

1410 

08 

1300 

OB 

3100 

&HR25 

EQU 

* 

• 

* 

CHARI 

9 

OB 

0010 

03 

1100 

OB 

1300 

OB 

1410 

OB 

1450 

D  3 

1360 

03 

1160 

OB 

1050 

DB 

1030 

OB 

1120 

DB 

1320 

OB 

3430 

•** 

• 

CHR26  EQU 

S 

;  CHARI  1 

OB 

0140 

09 

1130 

OB 

0120 

03 

3110 

CHR27 

EQU 

t 

• 

* 

CHAR  t 

• 

• 

08 

0240 

- 

03 

1230 

03 

0220 

03 

3100 

:•* 

5hR28 

EQU 

S 

• 

CHARI 

< 

03 

0460 

03 

1030 

03 

3400 

CHR29 

EQU 

5 

• 

• 

CHAR  I 

- 

03 

00  40 

03 

1440 

08 

3020 

03 

3420 

:*• 

3  HP  30 

EQU 

$ 

• 

• 

CHARI 

> 

03 

1430 

03 

3060 

CHR31 

EQU 

{ 

• 

♦ 

CHAR  I 

? 

OB 

0050 

03 

1160 

OB 

1260 

03 

1350 

OB 

1340 

DB 

1130 

OB 

1120 

DB 

0110 

□  3 

3100 

:♦* 

CHR  32 

EQU 

$ 

• 

• 

CHARI 

a 

C  3 

0330 

03 

1220 

03 

1130 

D  3 

1240 

OB 

1340 

OB 

1320 

OB 

1430 

OB 

1440 

DB 

1350 

OB 

1150 

DB 

1040 

03 

1020 

D  3 

1110 

03 

1310 

03 

3420 

:** 

ENQ1 

EQU 

5 

• 

***** 

***** 

ORG 

GCOOH 

:  ** 

:  ** 

PLOTS 

EQU 

s 

• 

it 

PUSH 

3 

PUSH 

0 

• 

PUSH 

H 

• 

• 

LXI 

B  *  C 

• 

t 

CON VI EM T 

CONSTANTS 

LHLD 

CURE  X 

• 

• 

* 

WITH  CJRFX  AND  RESETS  CUKR 

XCHG 

lhlo 

ENOx 

• 

• 

LOAuS  OE 

SHlD 

CJPPX 

• 

• 

« 

.-'N.. •  -  A  V\ 


n! 


M 


I 

■') 


• 

XCHG 

• 

CALL 

CMPAR 

:  COMPARES  ( 0£ )  TO  (HL) 

JC 

NEGDX 

:  endx <curkx 

• 

JN7 

POSOX 

:  ENO X>CURRX 
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University  of  California 
Los  Alamos  Scientific  Laboratory 
Attn:  Librarian  MS  362 

Sandia  Laboratories 

Attn:  DOC  CON  for  W.B.  Brown, 
Org.  1353 

Attn:  Tech  Library  Org.  3141 

Argonne  National  Laboratory 
Records  Control 

Attn:  DOC  CON  for  D.W.  Green 
Attn:  DOC  CON  for  LIR  SVCS  Rpts 
Sec. 

Attn:  DON  CON  for  G.T.  Reedy 

University  of  California 
Lawrence  Livermore  Laboratory 

Attn:  W.H.  Duewer,  L-262 
Attn:  J.  Chang,  L-71 

I'.S.  Energy  Rsch.  &  Dev.  Administration 
Division  of  Headquarters  Services, 
Library  Branch 

Attn:  DOC  CON  for  Class.  Tech. 
Library 


OTHER  GOVERNMENT 


US  ENERGY  RSCH,  AND  DEV.  ADMIN 

Division  of  Military  Application 
U.S.  Energy  Rsch.  &  Dev.  Admin. 
Attn:  DOC  CON 


Department  of  Transportation 
Office  of  the  Secretary 

Attn:  S.C.  Coroniti 

NASA 

Langley  Station 

Attn:  Tech.  Library 
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NASA 

Ames  Research  Center 

Attn:  N-245-3  R.  Whitten 

Department  of  the  Army 

Bal.  Missl.  Def.  Adv.  Tech.  Ctr. 

Attn:  W.O.  Davies 

Federal  Aviation  Administration 

Attn:  HAPP/AEQ -10 /James  W.  Rogers 

Central  Intelligence  Agency 

Attn:  ED/SI  RM  5G48  HQ  Bldg. 

Attn:  NED/OS  I-2G4R  HQS 

Department  of  Commerce 
National  Bureau  of  Standards 

Attn:  Sec.  Officer  for  M.  Krauss 
Attn:  Sec.  Off.  for  L.H.  Gevantman 

National  Oceanic  &  Atmospheric  Admin. 
Environmental  Research  Laboratories 
Department  of  Commerce 
Attn:  G.  Reid 
Attn:  E.  Ferguson 

Attn:  F.  Fehsenfeld 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


Science  Applications,  Inc. 


Attn: 

D.G.  Hopper 

Aero-Chem 

Research  Laboratories 

Attn: 

A.  Fontign 

Attn: 

H.  Pergament 

Aerodyne  Research,  Inc. 

Attn: 

F.  Bien 

Attn: 

M.  Camac 

Aerospace 

Corporation 

Attn: 

N.  Cohen 

Attn: 

H.  Mayer 

Attn: 

R.J.  McNeal 

Attn: 

T.D.  Taylor 

Attn: 

J.  Reinheimer 

Attn: 

R.D.  Rawcliffe 

Attn: 

R.  Herm 

Attn: 

P.  Kisliuk 

Battelle  Memorial  Institute 
Attn:  H.L.  LaMuth 
Attn:  STOIAC 

Brown  Engineering  Company,  Inc. 

Attn:  N.  Passino 

General  Research  Corporation 
Attn:  D.  Jones 
Attn:  J.  Ise,  Jr. 

California  At  Riverside,  University  of 
Attn:  J.N.  Pitts,  Jr. 

Attn:  A.M.  Winer 

California  At  San  Diego,  University  of 
Attn:  S.C.  Lin 

California  University  of  Berkley 

Attn:  Sec.  Off.  for  H.  Johnston 
Attn:  Sec.  Officer  for  Dept  of 
Chem. ,  H.L.  Strauss 

Calspan  Corporation 


Attn: 

C.E.  Treanor 

Attn: 

J.M.  Grace 

Attn: 

M.G.  Dunn 

Attn: 

W.  Wruster 

University  of  Colorado 

As  t  r o-Ge  ophys i cs 

Attn: 

J.B.  Pearce 

Colorado,  University  of 
Office  of  Contracts  and  Grants 

Attn:  G.M.  Lawrence,  LASP 

Concord  Sciences 

Attn:  E.A.  Sutton 

University  of  Denver 
Space  Science  Laboratory 
Attn:  B.  Van  Zyl 

University  of  Denver 
Denver  Research  Laboratory 

Attn:  Sec  Officer  for  D.  Murcray 
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AVCO-Everett  Research  Laboratory  Inc. 
Attn:  Tech.  Library 
Attn:  C.W.  Von  Rosenberg,  Jr. 

General  Electric  Company 
Space  Division 

Attn:  M.H.  Bortner,  Space  Sci.  Lab 

Attn:  J.  Burns 

Attn:  F.  Alyea 

Attn:  P.  Zavitsands 

Attn:  R.H.  Edsall 

Attn:  T.  Baurer 

Geophysical  Institute 
University  of  Alaska 

Attn:  J.S.  Wagner 
Attn:  Neal  Brown 

Lowell  University  of 
Center  for  Atmospheric  Research 
Attn:  G.T.  Best 


Lockheed  Missiles  and  Space  Company 
Attn:  J.  Kumer,  Dept  52-54 
Attn:  J.B.  Cladis,  Dept  52-12,  8202 
Attn:  B.M.  McCormac,  Dept  52-54 
Attn:  T.  James,  Dept  52-54 
Attn:  M.  Walt,  Dept  52-10 
Attn:  R.D.  Sears,  Dept  52-54 

Institute  for  Defense  Analysis 
Attn:  E.  Bauer 
Attn:  H.  Wolfhard 

Mission  Research  Corporation 
Attn:  D.  Archer 
Attn:  D.  Fischer 
Attn:  M.  Scheibe 
Attn:  D.  Sappenfield 
Attn:  D.  Sowle 

Photometries,  Inc. 

Attn:  I.L.  Kofsky 

Berkeley  Research  Associates 
Attn:  J.B.  Workman 

Physical  Dynamics,  Inc. 

Attn:  A.  Thomas 


General  Electric  Company 


Tempo-Center 

for  Advanced  Studies 

Attn: 

DASAIC 

Attn: 

W.S.  Knapp 

Attn: 

T.  St.eph.ens 

Attn: 

D.  Chandler 

Attn: 

V.R.  Strull 

Physics  International  Company 

Attn: 

DOC  CON  for  Tech  Library 

Pittsburg,  University  of  the  Comwlth 
System  of  Higher  Education 

Attn: 

W.L.  Fite 

Attn: 

M.A.  Biondi 

Attn: 

F.  Kaufman 

R&D  Associates 

Attn: 

R.  Latter 

Attn: 

R.G.  Lindgren 

Attn: 

B.  Gabbard 

Attn: 

R.  Lelevier 

Attn: 

A.L.  Latter 

Attn: 

F.  Gilmore 

Attn: 

H.J.  Mitchell 

Rand  Corporation 

Attn: 

C.  Crain 

Science  Applications,  Inc. 

Attn: 

D.A.  Hamlin 

Attn: 

D.  Sachs 

Attn: 

D.G.  Hopper 

Stanford  Research  Institute  International 

Attn. 

M.  Baron 

Attn: 

W.G.  Chesnut 

Technology  International  Corporation 

Attn: 

W.P.  Boqulst 

United  Technologies  Corporation 

Attn: 

R.H.  Bullis 

Utah  State  University 

Attn: 

D.  Baker 

Attn: 

K.  Baker 

Attn: 

C.  Wyatt 

Attn: 

A.  Steed 
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Physical  Science,  Inc.  Visidyne,  Inc. 


Attn: 

K.  Wray 

Attn: 

H.  Smith 

Attn: 

R. L.  Taylor 

Attn: 

J.W.  Carpenter 

Attn: 

G.  Caledonia 

Attn: 

T.C.  Degges 

Attn: 

C.  Humphrey 

Comma nder 

Attn: 

W.P.  Reidy 

Rome  Air  Development  Center 

Attn:  OSCA,  J.J.  Simons  Wayne  State  University 


Attn:  R.H.  Kummler 
Attn:  W.E.  Kaupplia 

Steward  Radiance  Laboratory 
Attn:  R.  Huppi 

Boston  College 

Space  Data  Analysis  Laboratory 
Attn:  E.R.  Hegblom 
Attn:  W.F.  Grieder 

Forrestial  Campus  Library 
Princeton  University 
Attn:  Librarian 
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